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ABSTRACT 


A Comparison of Experimental and Theoretical Results 
for Labyrinth Gas Seals. (May 1987) 

Joseph K. Scharrer, B.S., Northern Arizona University; 

M.S., Texas A&M University 

Chairman of Advisory Committee: Dr. Dara Childs 

The basic equations are derived for a two-control-volume model 
for compressible flow In a labyrinth seal. The flow Is assumed to be 
completely turbulent and Isoenergetic. The wall friction factors are 
determined using the Blaslus formula. Jet flow theory Is used for the 
calculation of the recirculation velocity In the cavity. Linearized 
zeroth and first-order perturbation equations are developed for small 
motion about a centered position by an expansion In the eccentricity 
ratio. The zeroth-order pressure distribution is found by satisfying 
the leakage equation. The circumferential velocity distribution is 
determined by satisfying the momentum equations. The first order 
equations are solved by a separation of variable solution. 
Integration of the resultant pressure distribution along and around 
the seal defines the reaction force developed by the seal and the 
corresponding dynamic coefficients. The results of this analysis are 
compared to experimental test results presented In this report. The 
results presented are for three teeth-on-rotor and three teeth-on- 
stator labyrinth seals with different radial clearances. The theory 
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compares well with the cross- coupled stiffness data Tor both seal 
types and with the direct damping data for a teeth-on-rotor labyrinth 
seal. For a teeth-on-stator labyrinth seal, the test results show a 
decrease In direct damping for an Increase In radial seal clearance, 
while the theory shows the opposite. 
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NOMENCLATURE 


Aj Cross sectional area of control volune (L a ); illustrated in 
figures (9) and (14) 

Bj Height of labyrinth seal strip (L); illustrated in figure (6) 

C Direct damping coefficient (Ft/L) 

C r Nominal radial clearance (L); illustrated in figure (6) 

D h Hydraulic diameter of cavity (L); introduced in equation (44) 

H Local radial clearance (L) 

K Direct stiffness coefficient (F/L) 

L Pitch of seal strips (L); illustrated in figure (6) 


NT Number of seal strips 
NC-NT-1 Number of cavities 
P Pressure (F/L a ) 

R Gas constant (L a /Tt a ) 

Rs, Radius of control volume 1 (L); illustrated in figure (6) 

Rs 2 Radius of control volune II (L); illustrated In figure (6) 

Rsu Surface velocity of rotor (L/t) 

T Temperature (T) 

Tp Tooth tip width (L); illustrated in figure (6) 

U, Average axial velocity for control volume I (L/t); Illustrated in 
figure (7) 

U 2 Average axial velocity for control volume II (L/t); illustrated in 
figure (7) 

W,i Average circumferential velocity for control volune I (L/t); 
illustrated in figure (7) 


W aI Average circunferentlal velocity for control volune II (L/t); 
Illustrated In figure (7) 

Wo l Average circunferentlal velocity In the Interface between control 
volunes I and II (L/t); Introduced In equation (-M0) 
a,b Radial seal displacement components due to elliptical whirl (L); 
Introduced in equation (73) 

ar Dimensionless length upon which shear stress acts on rotor; 
introduced in equation (H2) 

as Dimensionless length upon which shear stress acts on stator; 
introduced in equation (il2) 

c Cross coupled damping coefficient (Ft/L); in equation (ill ) 
e D Displacement of the seal rotor from centered position (L) 
k Cross coupled stiffness coefficient (F/L); in equation (ill) 
m Leakage mass flow rate per circumferential length (M/Lt) 

mr, nr, ms, ns Coefficients for friction factor; introduced in 
equation (I43) 
t Time (t) 

v Total velocity (L/t); introduced in equation (I18) 
to Shaft angular velocity (1/t) 

0 Shaft precesslonal velocity (1/t) 

P Density of fluid (M/L») 
v Kinematic viscosity (L*/t) 
e - e 0 /C r Eccentricity ratio 

e Turbulent viscosity (Ft/L* ); introduced in equation (9) 

* - 3. 1il1592 

Y Ratio of specific heats 


Subscripts 


0 Zeroth-order component 

1 First-order component, control volume I value 

2 Control volume II value 
1 1-th chamber value 

J Value along the dividing streamline 
x X-direction 
y Y-dlrection 
r Reservoir value 


s Sump value 


CHAPTER I 


INTRODUCTION 

The problems of Instability and synchronous response In 
turboraachlnes have arisen recently because of the trends in design 
toward greater efficiency with higher performance. To achieve these 
design goals, the machines are designed for higher speeds, larger 
loadings, and tighter clearances. In order to achieve the higher 
speeds, rotors frequently traverse several critical speeds (speeds 
which coincide with the rotor's damped natural frequency). The 
characteristics of synchronous reponse, when the rotor vibrates at a 


frequency coincident 

with the 

running 

speed. 

are 

such 

that 

the 

vibrational amplitude 

reaches a 

maximum 

at each 

critical 

speed. 

In 

order to limit the peak 

synchronous vibration levels, 

damping must 

be 


introduced into the rotor system. As loadings are increased and clear- 
ances decreased, fluid forces increase and can lead to unstable or 
self-excited vibrations. This motion is typically subsynchronous, 
which means that the rotor whirls at a frequency less than the rotating 
speed, and occurs with large amplitudes which grow as running speed 
increases. This situation can also be improved by adding damping to 
the rotor system, which would help curb the growth of the amplitudes. 
One of the rotordynamic force aechanlsms which plays a role in self- 
excited vibration and synchronous response is that of the forces devel- 
oped by labyrinth seals. 
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A limited amount of experimental data has been published to date 
on the determination of stiffness and damping coefficients for 
labyrinth gas seals. The first published results for stiffness 
coefficients were those of Wachter and Benckert [1,2,3]. They 
investigated the following three types of seals: a) teeth-on-stator, b) 
teeth on the rotor and stator, and c) teeth on the stator and steps or 
grooves on the rotor. These results were limited in that the pressure 
drop was small, much of the data were for nonrotating seals, no data 
were presented for seals with teeth on the rotor, the rotor speed was 
limited, and tests where rotation and inlet tangential velocity existed 
simultaneously were very scarce. 

The next investigation was carried out by Wright [4], who measured 
an equivalent radial and tangential stiffness for single-cavity seals 
with teeth on the stator. Although for a very limited and special 
case, Wright's results do give insight into the effect of pressure 
drop, convergence or divergence of the clearance, and forward or 
backward whirl of a seal. These results could be reduced to direct and 
cross-coupled stiffness and damping, hence, they are the first 
published damping coefficients for teeth-on-stator labyrinth seals. 

Brown and Leong [5] investigated the same seal configurations as 
Wachter and Benckert, in an effort to verify and extend their work. 
Their results Include variations of pressure, geometry, rotor speed, 
and inlet tangential velocity. Although the investigation was 
extensive, the published results are limited beca\»ee of the lack of 
Information concerning operating conditions for the various tests. 

Childs and Scharrer [6] Investigated geometrically similar teeth- 
on-rotor and teeth-on-stator labyrinth seals for stiffness and damping 


coefficients up to speeds of 8000 cpro. Kanemltsu and Ohsawa [7] 
Investigated multistage teeth-on-stator and interlocking labyrinth 
seals up to speeds of 2*400 cpm. They measured an effective radial and 
tangential stiffness while varying the whirl frequency of the rotor. 
These data could be reduced to stiffness and damping coefficients. 
Hlsa et al. [8] Investigated teeth-on-stator seals with 2-U teeth and a 
teeth-on-stator seal with steps on the rotor up to speeds of 6000 cpm. 
These data only Included static tests for direct and cross— coupled 
stiffness using steam. 

In the area of theoretical analysis of labyrinth seals, there is 
much more published information. The first steps toward the analysis 
of a labyrinth seal were taken by Alford [9] , who neglected 
circumferential flow and Spurk et al [10] who neglected rotation of the 
shaft. Vance and Murphy [11] extended the Alford analysis by 
introducing a more realistic assumption of choked flow. 

Kostyuk [12] performed the first comprehensive analysis, but 
failed to include the change in area due to eccentricity which is 
responsible for the relationship between cross-coupled forces and 
parallel rotor displacements. Iwatsubo [13. 1 1 *] refined the Kostyuk 
model by including the time dependency of area change but neglected the 
area derivative in the circumferential direction. Kurohashl [15] 
incorporated dependency of the flow coefficient on eccentricity into 
his analysis, but assumed that the circumferential velocity in each 
cavity was the same. Gans [16] improved on the Iwatsubo model by 
introducing the area derivative in the circumferential direction. 
Martinez-Sanchez et al. [17] produced results similar to Gans, but used 


empirical flow coefficients to improve their results. Childs and 
Scharrer [18] improved the Iwatsubo solution by using a modified set of 
reduced governing equations and an efficient solution technique. 
However, their solution continued to neglect axial velocities. 

Hauck [19] introduced the use of multi-control-volume analysis in 
the study of labyrinth seals by applying the equations of Impulse and 
"balance of moments" to a three-control-volume model. These equations 
were written in the axial direction only and neglected the effects of 
rotor speed. Fujikawa et al. [20] Introduced the use of two control 
volumes into the analysis of labyrinth seals, but their analysis, which 
neglects the axial velocity components, was heavily dependent on 
empirical information which is not customarily available. Finally, 
Jenny et al. [21] used the two-control-volume approach in conjunction 
with a two dimensional solution to the Navier-Stokes equations to 
account for the free shear stress between the Jet flow and the cavity 
flow. However, they neglected the recirculation velocity in the 
cavity, assumed the flow to be incompressible, and their free shear 
stress relation required a correction factor to fit the experimental 
data. Further, the present author obtains different signs in the 
expansion of the continuity equation and different perturbation 
equations. These discrepancies are explained in detail in the 
following section. 

The most extensive comparison of analytical predictions and 
experimental results was carried out by Scharrer [22] using the theory 
of Childs and Scharrer [18] and the results of Childs and Scharrer [6], 
This comparison showed that the theory [18] predicts cross-coupled 


stiffness reasonably well, but undcrpredicts direct stiffness, direct 
damping, and cross-coupled damping. 

In reviewing the state of the art In labyrinth seal 

experimentation and analysis, it becomes clear that there is a need for 
(a) an improved theory for the prediction of damping coefficients, (b) 
more extensive testing of teeth-on-rotor seals, and (c) test results 
showing the effects of change of radial seal clearance and higher 
speeds on stiffness and damping coefficients. This report will 

describe the revised test facility and program designed to measure the 
forces developed in a gas labyrinth seal. Some results, showing the 
effect of radial clearance change and higher rotor speed on teeth-on- 
rotor and teeth-on-stator labyrinth seals, will be presented. Also, a 
new analysis, which more accurately describes the physics of the flow 
field, will be presented and a comparison made between the theoretical 
and experimental mass flow rate, pressure distribution, and seal 
coefficients for the new results presented here. 

The major contribution of this report is a new analysis which 
incorporates the recirculating velocity in the cavity into the shear 
stress calculations. In addition, the analysis is based on a close 
comparison with the CFD results of Rhode [23,2*0. The CFD results 
were reinforcement for the assumptions and modelling of the shear 
stresses and velocity profiles. Stoff [25] carried out a comparison of 
CFD and experimental results for incompressible flow in a labyrinth 
seal. However, not enough data was provided in the paper to use the 
results for a comparison in this report. 


CHAPTER II 


THEORETICAL DEVELOPMENT 

SEAL ANALYSIS OVERVIEW 

As related to rotordynamlcs, seal analysis has the objective of 
determining the reaction forces acting on the rotor arising from shaft 
motion within the seal. There are two linearized seal models, expressed 
in terms of dynamic coefficients, which have been suggested for the 
force-motion relationship. For small motion about an eccentric 
position, as shown in figure 1, the relations of equation (1) have been 
proposed . 


Kxx(c 0 ) KxyiEoH^Xl rCxx(eo) C X y(e 0 )*UX| 

K yx(eo) Kyy(eo)J(Y/ LCy X (eo) Cyy(eo)J(Y» 

where the dynamic coefficients {Kxx ,Kyy ,Cxx ,Cyy) and {Kxy,Kyx, Cxy.Cyx} 
represent the direct stiffness and damping and the cross-coupled 
stiffness and damping, respectively. These coefficients are functions 
of the equilibrium eccentricity ratio, c 0 - e 0 /C r . where e 0 is the 
displacement of the rotor from the centered position and C r is the 
nominal radial clearance. The cross-coupling terms result when motion 
in one plane results in a reaction in a plane orthogonal to it. These 
cross-coupling terms depend on the magnitude and direction (with 
respect to the rotor*s rotation) of the fluid*s circumferential 
velocity. This velocity may exist at entry to the seal or may develop 
as the fluid passes through the seal. The cross-coupled stiffness term 



\\ 



Small motion of a' seal rotor about an eccentric position 
u) is the rotor 9 pin speed, /I is the precessional orbit 
frequency. 



Small motion of a 9eal rotor about a centered position, 
oj is the rotor spin speed , fl is the precessional orbit 
frequency. 



usually produces a destabilizing force component, and Is therefore of 
considerable Interest. The cross-coupled damping term is generally 
much less significant than the cross-coupled stiffness term with 
respect to stability. The second linearized seal model Is applicable 
for small motion about a centered position, as shown In figure 2. The 
form of the model is 


m : :]|:h: :k:i 


( 2 ) 


where the dynamic coefficient matrices are skew-symmetric. This model 
Is used in the analysis which follows. 

Preamble 

The flow in a labyrinth seal has been shown by experiment 
[1i»] and calculation [23] to be comprised of two flow regimes: a Jet 
flow region In the leakage path and a recirculating velocity region in 
the cavity itself(see figure 3). The first attempts at analysis of 
this system neglected the axial velocity components In the flow and 
concentrated on the circumferential components. This was the single 
control volume approach, used in refs [9~18]. In an attempt to Improve 
upon the results of these analyses the two-control- volume approach was 
Introduced, see refs [20,21]. These analyses Incorporated the axial 
velocity of the Jet flow into the solution but not the recirculating 
velocity component of the cavity flow. The results from Jenny et al. 
[21] showed substantial improvement in the prediction of stiffness and 
damping coefficients, but, in the end, correction factors had to be 
Incorporated into the calculation of the shear stress to improve the 


correlation with test data. 
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Fig. 3 Flew pattern in a labyrinth seal cavity. 
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Fig. 4 TVo-oontrol-voliine model with 
recirculation velocity, U 2 . 
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This report introduces the calculation of the recirculation 
velocity into the analysis. The model for the recirculation velocity, 
U,, used here is illustrated in figure A. This velocity component is 
Important in the calculation of the cavity shear stresses. The focus is 
on the shear stresses, because experimental results [6] have shown that 
the stiffness and damping coefficients are very sensitive to the 
circumferential velocity in the seal. In the control volume analysis 
to be presented, the solution to the circunferential momentum equation 
yields the circumferential velocity in the seal. An improvement in the 
shear stress calculation will yield an improvement in the calculation 
of the stiffness and damping coefficients. 

Before proceeding with the solution development, the approach 
taken in modelling the flow will be discussed. As mentioned 
previously, the flow in a labyrinth seal is known to have two distinct 
regions: a Jet flow region in the leakage path and a recirculating flow 
region in the cavity ltself(see figure 3). Therefore, a two-control- 
volume model seems appropriate. The choice is between the "box-in-a- 
box" model (see figure 5) of Jenny et al [21] or a more conventional 
model with a control volune for the Jet flow and one for the 
recirculating flow in the cavity, as shown in figure A. The two- 
separate-control- volume model was chosen, since It is suggested by the 
known physics of the flow. The flow enters the seal and separates into 
two distinct flow regions which are separated by the dividing 
streamline. 

The final question is whether the control volunes should be 
defined using a geometric boundary or using the dividing streamline as 
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Fig. 5 The "box-in-a-box" control volume 
model of ref.[2l] 


the boundary. The dividing streamline approach seems, at first, to be 
the obvious choice. The governing equations would be simplified by the 
restriction of no flow across a streamline, the free shear stress 
relations are derived for flow along the dividing streamline, and the 
solution for the velocity of the recirculating flow may be derived for 
flow along the dividing streamline. Despite these advantages, the 
dividing streamline approach was not used, however, it will now be 
reviewed and a method of solution discussed. The geometric boundary 
approach and solution is provided in the following section. 

DIVIDING STREAMLINE APPROACH 
Assumptions 

The following assumptions are used in deriving the governing 
equations: 

1) The fluid is considered to be an ideal gas. 

2) Pressure variations within a chamber are small compared to the 
pressure difference across a seal strip. 

3) The lowest frequency of acoustic resonance in the cavity is much 
higher than that of the rotor speed. 

A) The eccentricity of the rotor is small compared to the radial seal 
clearance. 

5) Although the shear stress is significant in the determination of 
the flow parameters (velocity etc.), the contribution of the shear 
stress to the forces on the rotor are negligible when compared to 
the pressure forces. 

6) The cavity flow is turbulent and lsoenergetic. 

7) The recirculation velocity, 1),, is unchanged by viscous stresses as 
it swirls within a cavity. 


Procedure 


The following analysis is developed for the teeth-on-rotor "see- 
through" labyrinth seal shown In figure 6. The continuity and 
circumferential momentum equations are derived for the two-control- 
volume model shown In figures 7 through 11. A procedure is discussed 
for determining the approximate location of the dividing streamline and 
the perturbation of the dividing streamline for small motion about a 
centered position. 

Continuity Equations 

Figures 7 and 8 show the control volumes defined by the dividing 
streamline. These control volumes have a unity circumferential width. 
Their continuity equations are: 

3pA , 9 pW , A , 


* + m i*1 - mi • 0 

3t Rs | 96 

3pA a 3pW 2 Aj 


(3) 


II: 

where the control 
are defined by 


♦ - o 

3t Rs 2 36 

volume areas, A, and A 2 , are shown in 


(*0 

figure 9 and 


L L 

A, » LCr + / ydx ; A 2 - LB - / ydx ( 5 ) 

0 0 


The following momentum equations for control volumes I and II are 
derived using figures 10 and 11 which show the pressure forces and 
shear stresses acting on the control volumes. It Is important to note 
that the recirculation velocity, U 2 , is Included in the shear stress 
definitions used In equations (6) and (7). These definitions are 
developed in a subsequent section of this chapter. 
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Fig. 10 Forces on control volumes 



Fig.ll Pressure forces on control volume I 
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Streamline Location 

The main difficulty in obtaining a solution to the above equations 
( 3 ) through ( 7 ) is in the determination of the location of the dividing 
streamline. The streamline definition, y(x), must be found to 
determine the control- volumes areas. A, and A a , defined in equations 5. 
There is no known solution for the location of a dividing streamline 
for the three dimensional flow field found in a labyrinth cavity. An 
approximation for the location of the dividing streamline can be 
obtained using the theory for the flow of a two-dimensional, turbulent, 
lsoenergetlc, half-infinite Jet. Figure 12 shows the model for this 
theory. The flow is assumed to enter with one velocity component, in 
the x-di recti on, and spread into the cavity, developing a y-coraponent 
of velocity. This model does not account for the circixnferential velo- 
city component, which is the same order of magnitude as the axial velo- 
city, in a labyrinth seal flowfield. The solution procedure involves 
solving the inf lnltessimal form of the x-momentum equation for the 
dimensionless velocity profile and then solving the integral form of 
the continuity and momentum equations for the location of the dividing 




streamline. A complete discussion of this theory can be found In Korst 
et al. [26]. 

The following Is a derivation of the equations necessary to 
determine the location of the Jet dividing streamline. The following 
derivation uses the assumption that the curvature in the dividing 
streamline Is small. The lnflnltessimal form of the continuity equation 
for the flow illustrated in figure 12 is: 

3(pu)/3x ♦ 3(pv )/3y - 0 (8) 

where the x and y velocity components, u and v, respectively, are time 
averaged. The lnflnltessimal x-moroentum equation, which has been 
reduced using equation (8), Is: 

pu 3u/3x ♦ pv 3u/3y - 3(ep3u/3y)/3y ( 9 ) 

where e Is the apparent (turbulent ) kinematic viscosity. Since the flow 
illustrated in figure 12 is a quasl-one-dimenslonal Jet flow where 
there is little or no Initial vertical velocity component, equation (9) 
can be linearized using the following perturbation method: 

u - u, ♦ V" 5 v - V" ; p - Pl ♦ p* do) 

where |v"|«|U,| and |u"|«|U,|. The resultant equation is: 

p,U, 3u"/3x ♦ p,V, 3u"/3y - ep, 3*u"/3*y (11) 

where e-e(x) and the second terra (V, term) is considered small. The 
final form of the equation is: 

U, 3u"/3x - e 9 a u n /By a (12) 

The following dimensionless variables are introduced: 

+ - U/U, - 1+uVU, 

C - y/6 

* ■ */6 (13) 

* 

C - / e d*/(U,«) 

0 

where 6 is the initial boundary layer thickness shown in figure 12. 


E 


Equation (12) becomes: 


3$/H - 3 *4>/3c 5 


with the Initial conditions: 


♦ • ♦(O.t) - 0 for -• < c < 0 

♦ - ♦(0,c) - ♦o (t ) for 0 < c < 1.0 

♦ ■ ♦(0,t) - 1.0 for 1.0< c < • 

and boundary conditions: 

♦ - ♦(5.—) - 0 for i > 0 

♦ - ♦(£,•) - 1.0 for E > 0 

The solution to equation (HO for the above Initial and boundary 
conditions is: 


♦ - 0.5[1~erf (n p ~n)] ♦ 1//T / ♦[(rrBJ/npJe dB (15 

n-n p 

where 

h p - i/(2/fc) ; n - cnp 
x 

erf(x) - 2 // V f exp(-$*)dfi 
0 

erf(-x) - -erf(x) 

The apparent viscosity, e, can be expressed in terms of the apparent 
viscosity far from the mixing region, e,: 


where 


e • e,f(^) 

f(<p) — > 1.0 as — > • 


According to Prandtl, this can be rewritten as: 


with 


e, - •cb(x)[U niax - Umin] 
3e/3y - 0 


For a half-infinite jet, equation (17) is: 



e, - «cb(x)U, (18) 

where b(x) Is the width of the mixing region. Assuning that the mixing 


Increases linearly, l.e. 


b(x) - cx - ci|>6 


where c is a constant, equation (16) becomes: 


e » c4>6U,f(ip) 


( 20 ) 


Substituting this into £, from equation ( 1 3 ) yields: 

* 

£ ■ c / i|»r((|>) d(p (21 ) 

0 

Looking at a limiting case of equation (21): 
as x/6 — > • 

then — > • 

i "> - 

n p - 1/2 vT --> 0 

This limiting case Is for either no initial boundary layer, which Is a 

good assumption for labyrinth seals, or fully developed velocity 

profiles. Since np — > 0, the variable n Is now undefined. Llepman 
and Laufer [27] have defined n for this limiting condition using the 
following development. By definition: 

as 1 > — > - then f(4») — > 1.0 
Inserting this into equation (21) yields: 

£ - c**/2 


By definition: 


n “ cn p - ( ;/ 2 /T - c/(¥>/2 c) - y/ (x/2c ) 

Letting c - 1/(2 o*), yields the desired result: 

n - oy/x (22) 

where c Is the Jet spreading parameter. Korst and Tripp [28] used 
experimental data to find the following relation for 0: 

0 - 12.0 ♦ 2.758M, (for air) (23) 

Coertler [29] has shown that the dimensionless velocity, t, follows 

directly from equation (15) when n p — > 0: 

♦ - 0.5(1*erf (n)) ( 2 U) 

Equations (2iJ) is a solution for the dimensionless velocity profile, 
at any dimensionless position, n. The goal of this development is to 


determine the dimensionless dividing streamline position, nj . nj can 
be obtained by solving the Integral form of the continuity and x- 
momentum equations for the system shown in figure 12. 

Control Volume Analysis 

The coordinate systems and definition of the control surface are 
shown In figure 12. The (x,y) coordinate system Is the Intrinsic coor- 
dinate system while the (X,Y) coordinate system is the reference 
system. Equations (22) and (2^) are approximate relations; exact 
relationships. If known, would provide for conservation of momentum for 
the constant pressure mixing region. The reference coordinate system 
is the coordinate system in which momentum is conserved. The intrinsic 

coordinate system Is located with respect to the reference coordinate 

system by a control volune analysis utilizing the conservation of 
momentum principle for this constant pressure mixing region. The 
relationship between the coordinate systems normal to the Jet is: 

ym(x) - y-Y 

wlth y m (0) - 0. 

X-Momentum Equation 

The steady flow x-momentura equation for the Jet flow shown in 
figure 12, written for the reference coordinate system and expressed in 
the previously defined dimensionless variables is: 

R R 

/ pU* dY | - / pU 1 dY | (25) 

0 X-0 — X-X 

For the momentum equation, the lower control surface is located at 
This equation contains no surface forces. This is realistic for a 
labyrinth seal if location R is far from of the stator wall. Rewriting 
equation (25) for the intrinsic coordinate system: 
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/ pU a dy | ♦ / pU* dy | - ? p0 a dy j 

0 x-0 6 x-0 0 


(26) 


x-x 


Introducing the previously defined dimensionless variables, equations 
(13) and (22), equation (26) becomes: 

1,0 * n R-Ttn 

n P / ( p/p,)^ e dc ♦ nR - np - / (p/p,U* dn (27) 

0 — » 

Distance R is chosen such that: 

1 ~ ♦( dr ) <« 1.0 

Equation (27) becomes: 

1 • 0 a hR 

n P / (p/pi)*o dc ♦ hr - np - / (p/Pi)4>* dn + Tta (28) 

0 — • 

Applying the condition of no initial boundary condition (np > 0) 

equation (28) is: 


hR 


% - nR - / (p/PiH* dn 


(29) 


Continuity Equation 

The steady flow continuity equation, written for the reference 
coordinate system, is: 

R R 

(30) 


R R 

/ pU dY J - / pU dY | 

0 X-0 yj- ym x-X 


For the continuity equation, the lower control surface Is coincident 
with the Jet dividing streamline. Rewriting equation (30 ) for the 
Intrinsic coordinate system: 

.6 R R+ 1 


6 R R+ y 

/ pU dy | + / pU dy j - / plTdy | 

0 x-0 j x-0 yj x . x 


(3D 


Introducing the previously defined dimensionless coordinates and 
multiplying equation (3D by n p /6: 
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n p / (p/PlHo d c ♦ HR - np - / ( p/p 1 )♦ dn 4 r»n (3?) 

0 

Substituting the results of the momentum equation, equation ( 29 ), into 
equation (32) yields: 

n R 1.0 hr 

/ (p/p , )♦ dn - / (p/p, ) ( 1 — ♦) o dc ♦ / (p/p»)i* dn (33) 

nj 0 

Making the assumption of no Initial boundary layer (n p — >0), equation 
(33) becomes: 

n R nR 

/ (p/p,)$ dn - / (p/p»)4>* dn (34) 

nj 

The density ratio, (p/p,), for isoenergetic flow (constant temperature) 
is given as: 

p/p, - (1-Ca*)/(1-Ca*4> a ) (35) 

The final form of the continuity equation becomes: 

^R nR 

/ ($/[1-Ca*$* J ) dn - / (♦*/[1-Ca*$*3) dn (36) 

nj 

Where Ca is the Crocco number. The Crocco number is defined as: 

Ca* - (Y-1 )M*/ (2+(Y-1 )M*) (37) 

The Crocco number is a dimensionless velocity similar to the Mach 
number. The Crocco number uses the maximum lsentroplc speed of a gas 
while the Mach number uses the local speed of sound. The Mach number 
varies between 0 and • while the Crocco number has a range of 0 to 1. 

The solution to equation (36), the location of the dividing 
streamline, can be obtained by the following steps: 

0) Calculate the Mach number using the zeroth-order leakage value. 
The zeroth-order leakage is discussed in the next section. 

1) Calculate the Crocco number using equation (37). 
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2) Substitute equation (2*0 Into equation (36) and Integrate the 
error function. The value of the error function at the limits R and -- 
Is 1.0, leaving an equation In nj only. This Is solved for n j , which 
Is the dimensionless location of the dividing streamline. 

3) Use the straight line approximation, equation (22), to find y 
as a function of x. 

**) Insert y(x) from step 3 into equation (5) and calculate the 
areas of the control volunes. 

The above procedure yields the zeroth-order (centered) value for 
the areas. The problem is to find the values for a perturbation analy- 
sis. Since equation (2*1) is an error function, an explicit equation 


for the change In area for a perturbation in cleai'ance cannot be 
obtained. However, the above procedure could be carried out for a 
range of clearances in the neighborhood of the nominal clearance and an 
approximation for the change in area and a final result could be 
obtained . 

As noted at the beginning of this discussion, the advantages of 
the dividing streamline approach are that the free shear stress and 
recirculating velocity equations may be derived along the dividing 
streamline, and the governing equations are simplified by the condition 
of no mass flow across a streamline. The above solution procedure 
yields only an approximation for the location of the dividing 
streamline for a simplified (two-dimensional) flow while increasing the 
difficulty in obtaining a solution. Therefore, the advantages of the 
dividing streamline approach are outweighed by the difficulty in 
obtaining a solution. The geometric boundary approach and a 


complete solution will now be presented. 


GEOMETRIC BOUNDARY APPROACH 


Procedure 

The analysis presented here Is developed for the teeth-on-rotor 
"see-through" labyrinth seal shown In figure 5. The equivalent 
equations for the teeth-on-stator labyrinth seal are given in Appendix 
A. The continuity and circumferential momentum equations will be 
derived for the two-control- volume model shown in figures 13,14,15, and 
16. A leakage model will be employed to account for the axial flow. 
The governing equations are linearized using perturbation analysis for 
small motion about a centered position. The zeroth-order continuity and 
momentum equations will be solved to determine the steady state 
pressure, axial and circumferential velocity for each cavity. The 
first-order continuity and momentum equations will be reduced to 
linearly independent, algebraic equations by assuming an elliptical 
orbit for the shaft and a corresponding harmonic response for the 
pressure and velocity perturbations. The force coefficients for the 
seal are found by integration of the first-order pressure perturbation 
along and around the shaft. 

Continuity Equations 

The control volumes of figures 13 and 14 have a unity 


circumferential 

width. 

3pA, 

Their continuity equations are: 
3p w iAj . .. .. 

Is 

at 

♦ ♦ m i + 1 <- mi ♦ m r - 0 

Rs,3e 

II: 


dpA a 3pW,A 2 . 


* - m r . o 

3t Rs,d6 


For the teeth-on-rotor case, A,-LCr, A,-LB, Rs,-Rs, and Rs,-Rs+B. 


Fig. 13 Control Volumes with 
geometric boundary. 



Fig. 14 Isometric view of Control Volumes 
with geometric boundary. 





28 


Momentum Equations 

The following momentum equations for control volumes I and II are 
derived using figures 15 and 16 which show the pressure forces and 
shear stresses acting on the control volumes. 
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where ar and as are the dimensionless length upon which the shear 
stresses act and are defined for the teeth-on-rotor labyrinth by 

asi - i ari - (2Bi «• Li) /Li . (42) 

W D is the circumferential velocity between the control volumes. 

Various models for the stator wall shear stress were evaluated by 
comparison to CFD results of Rhode [23]. For a teeth-on-rotor labyrinth 
seal, the optimum model for the stator shear stress (rotor shear stress 
for a teeth- on- st at or seal) was obtained by using the equation of 
Glauert [30] for wall shear stress of a plane Jet Issuing fom a slot. 
However, this relation requires knowledge of the maximum axial velocity 
and its displacement from the wall. This information is not available 
In a control volume analysis. The next best model, by comparison to 
[23], is Colebrook’s formula [31], but this equation is not explicit in 
the friction factor and cannot be perturbed. Experience [32] has shown 
that the perturbation of the friction factor Is important in stiffness 
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Fig. 15 Forces on control volumes 
with geometric boundary. 



Fig. 16 Pressure forces on control volume I 
of geometric boundary model. 


calculations. The next best shear stress model Is based on the assump- 
tion that the shear stresses (for rotor and stator surfaces) are 
similar to those found in the pipe analysis of Blaslus C 333 . Blasius 
determined that the shear stresses for turbulent flow in a smooth pipe 
could be written as 

> mo 

t - J. pU m no /UrnD h) 

2 V v / 

where Um is the mean flow velocity relative to the surface upon which 
the shear stress is acting. The constants mo and no can be empirically 
determined for a given surface from pressure flow experiments. 
However, for smooth surfaces the coefficients given by Yamada [3 1 *] for 
turbulent flow between annular surfaces are: 

mo - -0.25 no - 0.079 

Applying Blasius' equation to the labyrinth rotor surfaces yields the 
following definitions for the rotor shear stress in the circumferential 
direction. Note that the recirculation velocity, U a , is included in 
the definition of the total velocity acting on the rotor. 

1 /'(Rs a w-W a ) i *U a a Dhi\ rar 

■ — p/ 1 (Rs a w-W a )'*U a i (Rs a w-W a ) nn I (^3) 

2 \ t 

where Dh a j i 8 the hydraulic diameter of control volume II, defined by 

Dh al - 2BL/(B*L) (U*0 

Similarly, the stator shear stress in the circumferential direction is: 

1 /✓W 1 «*U 1 J Dhi\ m3 

*s - j W, ns^ — j (H5) 

where Dh,j i 8 the hydraulic diameter of control volume I, defined by 

Dh»i . 2CrL/ (Cr+L ) (46) 

and the axial velocity U, is 

(47) 


U, - ro/pCr 


Figure 17 shows a comparison of the predictions from equation ( ) 
and CFD results for stator wall shear stress for seal A of table 1 . 
The recirculation velocity, U lf is undefined at this point. It will be 
discussed In the following section. Table 1 shows the seal geometries 
calculated by Rhode [233* The figure shows that the comparison Is very 
good. Similar results are obtained for the other seals of table 1. 
Figure 18 shows a comparison of rotor wall-shear-stress predictions 
from equation (^3)* CFD, and averaged CFD results for rotor wall shear 
stress for seal A of table 1 . The averaged CFD result is used here for 
comparison since the bulk flow model yields a single averaged result 
for cavity shear stress and is not capable of modelling the complex 
11 owiiSiu. ihe figure shows that the prediction of equation ( 1 * 3 ) is 
close to the CFD results. The dips in the CFD results are the lower 
corners of the cavity. Similar results are obtained for the other 
seals of table 1 . 

The flow across a labyrinth tooth is very similar to the flow of a 

turbulent Jet issuing from a slot. The problem with using Jet-flow 

results for labyrinth seals is that current Jet-flow theory only 

considers the flow of a Jet with a coflowing stream or a crossflowing 

stream, not both. In the following derivation, the relations given by 

Abramovich [ 35 ] for the velocity profile of a semi-contained, one- 

dimensional, turbulent Jet with a coflowing stream are assumed to apply 

for the two-dimensional labyrinth seal flow. According to Abramovich 

[353, the velocity profile for such a flow can be shown to fit the 

following function when compared to experimental results: 

1.5 2 
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CFD Results 


Fig. 17 A comparison of Theoretical and 

CFD results for stator wall shear 
stress. See table 1 for seal 
geometry . 
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CFD Avg. 

CFD Results 

► Theory 


X/ (L+2B) 

Fig. 18 A comparison of Theoretical and 
CFD results for rotor wall shear 
stress. See table 1 for seal 
geometry. 





Fig. 19 Model of semi-contained turbulent jet 


Table 1 . Seal geometries calculated by Rhode. 


A 

Seal 

B 

C 

0 

72.05*1 3mm 

72.05*1 3mm 

72. 05*1 3mm 

*» 1.780mm 

3.175mm 

3.175mm 

3.175mm 

0.889mm 

3.175mm 

3.175mm 

3.969mm 

0.8585mm 

0.35mm 

0.35mm 

0.35mm 

0.15mm 

0.406*lmm 

O.508mm 

0.508mm 

0.2159mm 



where the coordinate y, the mixing thickness b, and the boundary layer 
thickness y, are defined in figure 19. The relationship between the 
boundary layer thickness and the mixing thickness was found [35] by 
comparison to experiment to be: 

y,/b - 0.584 - 0.13*»(v a /v,) (49) 

Once the velocity ratio across the dividing streamline, v a /v,, is 
found, equation (49) reduces to a constant. The total free shear stress 
is found using Prandtl's mixing length hypothesis [36]: 

T Jt - pi* | 3v/3y | Ov/3y) (50) 

where the nixing length, 1, for a labyrinth seal, has been determined 
from the calculations of Rhode [23] to be: 

l - 0.275b (51 ) 

Table 1 shows the seal geometries calculated by Rhode [22], The mixing 
length, l, given in equation (51) is the most sensitive factor in this 
solution. The large magnitude of the mixing length shows the high 
turbulence level of the labyrinth flow as compared to similar flows. 
The typical values given for the mixing lengths of rectangular and 
round Jet flows, in one dimension, are in the range of 0.07 to 0.09. 
Without the CFD results, one of these values would have to be used and 
the results of using l in the range [0.07,0.09] would have been disap- 
pointing. 

Jenny et al. [21] used a 2-D CFD code to obtain a correlation for 
1/b as a function of clearance and tooth geometry. Their relation is 
shown below for the teeth-on-rotor case: 

l/b - 0.055(1*1 .03Cr/L+0.08/Hs7E) (52) 

However, their shear stress relation neglected the recirculating 
velocity component, U,. Upon comparison with the data of Rhode [23], 


the mixing length ratio, t/b, was found to be relatively constant when 
the shear stress is calculated using all velocity components. 

Substituting the differentiated version of equation (i»8) and equa- 


tion (51) Into equation (50) yields an expression for the total free 
shear stress. At the Interface of the two control volumes (y*0), the 
total free shear stress Is: 

T Jt - 0.68 p|v a -v,|(v 1 -v I )[1-(y 1 /b]* 5 ] 2 (y 1 /b) (53) 
The circumferential component of the free shear stress is: 

T J - 0.68 p/(W a -W,W(U a -U,)’ (W a -W,)[1-(y a /b)* 5 ] 2 (y a /b) (5*0 
The circumferential component of the velocity at the Interface, W oi ^ la 
obtained from equation (*I8). 


W oi . w. ♦ (W a -W,)[1-(v./bl' 5 1 2 fv./h) rqm 
Equations (53*5*1,55) are all valid along the dividing streamline. 
Since the control volumes are defined geometrically and not by the 
dividing streamline, the shear stress calculated using the above 
equations is assumed to be close to that existing along the geometric 
boundary line. This is a good assumption considering that the angle of 
the dividing streamline from the horizontal has been found 
experimentally to be on the order of 6 degrees by several investigators 
[37.38]. 


The analysis to this point is incomplete in that the recirculation 
velocity, U a , and the relationship between the mixing thickness and the 
boundary layer thickness, y t /b, are undefined. In order to determine 
the recirculation velocity, U a , and subsequently y a /b, the analysis 
presented in the previous section dealing with the DIVIDING STREAMLINE 
APPROACH is used. Again, this analysis is valid along the dividing 
streamline, but is considered close enough to the values along the 


geometric boundary line. The final form of the continuity equation, 
equation (36), is rewritten here: 

n R HR 

/ (♦/[1-CaV]) dn - / (♦*/[1-CaVJ) dn (56) 

nj 

where Ca Is the Crocco number, n is the dimensionless coordinate, and 4 
Is the velocity ratio U/U,. The solution to this equation Is obtained 
by substituting equation (24) into equation (56) and solving for the 
dividing streamline coordinate, nj f for a given Crocco number. This is 
then inserted back into equation (24) and a value of i s obtained. 
The results of this solution procedure are tabulated in table 2, for 
air. For air (T-1.4) flowing in a labyrinth seal, the maximum possible 
Mach number is 1.0. Therefore, the maximum possible Crocco number is 
0.408 or Ca*-0.167. The range of solutions is: 

0.61632 < 4>j < 0.6263 

Using an average solution of - 0.62 gives a maximum error of less 
than ± 1%. The recirculation velocity at the interface is: 

u aj - 0.62U, (57) 

The only remaining problem is the numerical definition of y a /b. 
Looking back, equations (48) and (24) both describe the axial velocity 
profile in the Jet flowfleld. If the following observation is made 

v i/*i - 

then equation (57) can be substituted back into equation (49) yielding 
the following numerical definition for y t /b: 

y,/b - 0.584-0. 134 *j . o.50 

It is interesting to note that Jenny et al. [21] assumed that y a /b-0.5. 

Figure 20 shows a plot of the dimensionless axial velocity profile 
in the recirculation region for seal A of table 1 as calculated by 





Table 2. Tabulated solution to equation (56). 


Ca 7 


Ca 2 

h 

0.00000 

0.61632 

0.6R000 

0.67553 

0.05000 

0.61915 

0.72000 

0.68188 

0.10000 

0.62211 

0.76000 

0.68903 

0.15000 

0.62523 

0.80000 

0.69724 

0.20000 

0.62848 

0.84000 

0.70689 

0.24000 

0.63129 

0.86490 

0.713944 

0.28000 

0.63405 

0.88360 

0.719944 

0.32000 

0.63725 

0.90250 

0.726834 

0.36000 

0.64047 

0.92160 

0.734949 

0.40000 

0.64387 

0.94090 

0.744883 

0.44000 

0.64748 

0.96040 

0.757869 

0.48000 

0.65132 

0.98010 

0.777432 

0.52000 

0.65543 

0.992016 

0.798766 

0.56000 

0.65979 

0.998001 

0.823427 

0.60000 

0.66462 

1.000000 

1.000000 

0.64000 

0.66982 




Rhode [23]. This profile Is for the center of the recirculation region 
to the top of the labyrinth tooth. The intersection of the two dashed 
lines is the location and value of the theoretical recirculation velo- 
city as calculated using equation (57) and the assumption that the 
dividing streamline makes an angle of 6 # with the horizontal. The 
agreement is excellent. Equation (57) is actually the velocity at the 
Interface of the two control volumes. The velocity components used in 
the shear stress equations are all average velocity components. To be 
consistent, the average recirculation velocity must be used. The CFD 
results show that the velocity distribution is parabolic in nature. 
Integrating this yields: 


if _ r\ nnCi l 
v j v • C.VU v } 


reft 

W'-' / 


Reduced Equations 

The solution of the governing equations can be simplified by 
reducing the number of equations by one. This reduction is accomplished 

by using equation (39) to eliminate m r from the other equations. The 
continuity equation for control volume I becomes: 

3pA, 3pW,A, . . 9pA a 3pW a A a 

* ♦ n i-H - ml ♦ ♦ - 0 (59) 

3t Rs ,3e 3t Rs a 30 

If equation (59) times the circumferential velocity, W,, is now 

subtracted from equation (40), the following reduced form of the 

momentum equation for control volume I is obtained: 


3W, pW ,A | 3W, [3pA a 3W a A a pl 

pA, ♦ ♦ — ♦ < w oi-W,j 

3t Rs, 30 |_ 3t Rs a 30j 

A, 3P1 

4 n l(W,i-W,l-l ) ♦ xjiLl “ xsiasiLi 

Rs, 30 


( 60 ) 


Similarly, if equation (39) times the circumferential velocity, W a , 


is 


AO 


subtracted from equation (*U), the reduced momentum equation for 
control volume II Is obtained. 

aw, pWjA , aw, rapA, apw,A,*| 

P*a““ 4 ' 4 j 4 I ^ii~Woi ) (6 1 ) 

at rs, ae j_ at Rs,aeJ 

a, apj 

" 4 t JlLl - TrlariLl 

Rs, 30 

The number of variables is reduced by using the Ideal gas law to 
eliminate the density terms. 

Pi - plRT (62) 

This concludes the development of the governing equations for the new 
analysis presented in this report. The following is a discussion of 
the analysis of Jenny et al. [21]. 

The theory of Jenny et al. [21] 

The theory of Jenny et al. [21] has shown consistently good 
agreement with measured test results [39] in predictions of cross- 
coupled stiffness and direct damping. The author had hoped to program 
their solution and make direct comparison to the present theory; 
however, as outlined below, unresolvable difficulties arose in deriving 
the published equations of [21]. 

The theory of Jenny et al. [21] was derived for the "box-ln-a-box" 
control volume configuration Illustrated In figure 5. Thus, a direct 
comparison of their equations with those presented In this report is 
not feasible. However, a review of the development of their governing 
equations Is of Interest. 

The following convention will be used for the control volumes in 
figure 5: the large control volume Is control volume I and the small 
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control volume is control volume II. The continuity equations for the 
control volumes shown in figure 5 are: 

Continuity I 

apH,*, 3pW,A, 3p(A,*A a ) 

" 4 - mi - 0 (J1) 


Continuity II 


Rsde 


apWjAj 


Rsde 


at 


apA, 


“ ®ri - 0 


(J2) 


Rsae at 

The following assumptions are used by Jenny et al. [21] to simplify 
equations (J1) and (J2): 

a) the flow is incompressible (p - constant), 

b) Aj + i - bi t an( j 

c) the area of the control volume II Is constant. 

The first assumption seems questionable, since this is a compressible 

flow solution, and quite often the flow In a labyrinth seal achieves 

Hach 1 at the exit. Assumption (b) Is a valid assumption for the 

zeroth-order, steady flow solution, but It Is questionable for the 

first-order, unsteady flow solution for an orbiting rotor. Using the 

chain rule for the expansion of partial derivatives and the above 

assumptions, equations (J1) and (J2) become: 

Continuity I 

3W, 3W, acr 3(A,*A a ) 

Aj ♦ A, — ♦ W,L ♦ Rs - 0 (J3) 

ae ae ae at 


Continuity II 


aw t 

pA,'— - Rsm r i . o 

ae 


(JH) 


The equations given by Jenny et al. [21] are: 

Continuity I 

3W a aw, aCr 3(A,*A a ) 

A a — ♦ A, W,L - Rs - 0 

ae ae ae at 


(J5) 


Continuity II 


SW a 

pA, n r i - 0 (J6 ) 

de 

The difference between equations (J3) and (J5) la In the sign of the 
third and fourth terms. The second and third terms in equations (J3) 
and (J5) originate from the same partial derivative, but have opposite 
signs. This author could not arrive at the same conclusion using the 
chain rule. The difference between equations (J4) and (J6) is the 
radius, Rs, in the second term. This may or may not be a problem since 
the radial mass flow term, m r i t j 8 no t defined by Jenny et al. [21]. 

The author agreed with the derivation of the. momentum equations 
for the control volumes shown in figure 5 except for the aforementioned 
assumptions and the following discrepancies: 

(a) the axial velocity component is incorporated into the defini- 
tion of the stator wall shear stress, but neglected in the definition 
of the Reynold's number which is used to calculate the friction factor 
term in the shear stress relation. 

(b) the perturbation of the friction factor is Ignored. This term 
has been shown [32] to be Important in the solution for rotordynamlc 
coefficients. 

(c) the leakage equation is a global leakage equation. This means 
that local perturbations for a cavity can not be found from this equa- 
tion. Jenny et al. [21] perturb this global equation for clearance. 

(d) the carryover coefficient definition used in the leakage 
equation is a global equation and cannot be perturbed. 

(e) the flow coefficient used in the leakage equation was obtained 
from a plot of empirical data. No explanation was given for the method 


used to obtain the dervatives of the flow coefficient used in the per- 
turbation equations. 

The aforementioned problems prevented the author from obtaining a 
solution based on the theory of Jenny et al. [21]. Regrettably, no 
direct comparison between it and the theory presented In this paper was 
possible. This completes the discussion of the theory of Jenny et al. 
[21]. The following is a discussion of the solution procedure for the 
new analysis presented in this report. 

Leakage Equation 

To account for the leakage mass flow rate in the continuity and 
momentum equations, the following model was chosen. 


m l 


P»1 Pa Hi 


2 2 
p l-l ~ Pi 
RT 


( 63 ) 


where the kinetic energy carryover coefficient is defined by 

Vermes [*»()] for straight through seals as: 

1/2 

Ua * 1/[1~a] (6*1) 

where 


a - 8.52/((Li~T P i)/Cr*7.23) 

and is unity, by definition, for the first tooth of any seal and all 
the teeth in interlocking and combination groove seals. This 

definition of the carryover coefficient is a local coefficient which 
can be perturbed in the clearance. The previous analyses by Childs and 
Scharrer [6] and Jenny et al* [21] used a global definition which could 
not be perturbed. 

The flow coefficient is defined by Chaplygin [*I1] as: 


v i i 


v 


2 

«*2-5s 1 *2si 



where, si 


-1 


( 65 ) 



This flow coefficient yields a different value for each tooth along the 
seal as has been shown to be the case by Egli [42], For choked flow, 
Fliegner's formula C ^33 will be used for the last seal strip. It is of 
the form: 

. 0.510m, 

“NC - — — PNC HnC (66) 

✓ET 

PERTURBATION ANALYSIS 

For cavity 1, the continuity equation (59), momentum equations 
(60,61) and leakage equation (63) are the governing equations for the 
variables W a i, Pi, mj* A perturbation analysis of these equations 

is to be developed with the eccentricity ratio, e - e Q /Cr, selected to 
be the perturbation parameter. The governing equations are expanded in 
the perturbation variables 

P 1 - Pol ♦ e P ti Hi - Cri ♦ e H, 

W ii - W ,oi ♦ e W .,1 Ai - Ao ♦ e LH, 

W *i “ H ,oi ♦ e W 2I i 

where e ■ e Q /Cr Is the eccentricity ratio. The zeroth-order equations 
define the leakage mass flow rate and the circumferential velocity 

distribution for a centered position. The first-order equations define 

the perturbations In pressure and circumferential velocity due to a 
radial position perturbation of the rotor. Strictly speaking, results 
of a first order analysis are only valid for small motion about a 
centered position. 

Zeroth-Qrder Solution 

The zeroth-order leakage equation is 

*1*1 " *1 “ (67) 

and is used to determine both the leakage-rate m 0 and pressure 
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distribution for a centered position. The leakage-rate and cavity 
pressures are determined iteratively, in the following manner. First, 
determine whether the flow is choked or not by assuming that the Mach 
number at the last tooth is one. Then, knowing the pressure ratio for 
flow at sonic conditions, the pressure In the last cavity is found. 
The mass flow can be calculated using equation (66). Working 
backwards towards the first tooth, the rest of the pressures can be 
found using equation (63). The final pressure calculation will result 
in the reservoir pressure necessary to produce the sonic condition at 
the last tooth. If the actual reservoir pressure is. less than this 
value, then the flow is unchoked. Otherwise, it is choked. If the 

flow is choked, a similar procedure is followed, but now the pressure 
in the last cavity is guessed and a mass flow rate calculated using 
equation (66). The remaining pressures are calculated using equation 
(63). This is repeated until the calculated reservoir pressure equals 
the actual reservoir pressure. If the flow is unchoked, the pressure 
in the first cavity is guessed and a mass flow rate calculated using 
equation (63). The remaining pressures are calculated with the same 
equation. This procedure is repeated until the calculated sump 
pressure equals the actual sump pressure. 

The zeroth-order circumferential-momentum equations are 

m o(Wjoi-W,oi-1 ) • (Tjio~Tsio asi)Li (68) 

T JoiLi - Trol ariLi (69) 

From calculated pressures, the densities can be calculated at each 
cavity from equation (62), and the only unknowns remaining in equations 
(68) and (69) are the circumferential velocities W io j and W,oi. Given 
an inlet tangential velocity, a Newton-root-f indlng approach can be 
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used to solve equations (68) and (69) for the 1-th velocities, one 
cavity at a time; starting at the first cavity and working downstream. 
First-Order Solution 

The governing first-order equations (70,71,72), define the 
pressure and velocity fluctuations resulting from the seal clearance 
function. The continuity and momentum equations follow In order: 


8P, t 

0,1 — • 0.1 


3P,i 


3W , ,1 


♦ C , 1 - 

36 * 36 


♦ G„1 


3W,,i 

36 

3H.J 


+ G,i P,i 

3H,i 


(70) 


♦ G .i P,i-1 ♦ G,i P,i*i - - G.i H,l - G,1 - G le l ^ 


X U 


Will 


3t 

x »lPol 3W al l 
Rs 2 30 


*.i 


3W 2 ,i 
3t 


3P.1 


,iw l0 l 3W , , i r X # iW a of| 3P,i 3P,i 

— — * l_ x -‘ * — J jr * x -‘ ir ,7 ’ 

+ X Hl Pll ♦ X S 1 P,l-1 ♦ X.i w tl l ♦ X 7 1 w al l 
"® w »:i-1 " X.i H ,1 

pfiiPoi TUW.oilSWj.i r Y a iW a oi] 3Pjl 

♦ ♦ Y,1 ♦ (72 

|_ Rs a Rs a J 36 |_ Rs a J 36 


^ + Y h 1 P,i ♦ Y.l H a ,1 ♦ Y.l P,i-l ♦ Y 7 1 W.,1 - Y.i H,i 
where the Xj* 3 t Yi's, and Gi's are defined In Appendix B. These 
perturbation equations are very different from those of Jenny et ah 
[21], because their analysis neglects pressure perturbations in the 
leakage and shear stress equations and assumes that the density is 
constant. 

If the shaft center moves in an elliptical orbit, then the seal 
clearance function can be defined as: 

cH, - -a coswt cose “b slnwt sine 


- -a [cos (8-wt) ♦ cos (e*u>t)] - b [cos (0-wt) - cos(6+wt)] 
2 2 


(73) 


The pressure and velocity fluctuations can now be stated in the 
associated solution format: 
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4 4 - _ 

p >i “ Pclcos(e*wt) * P S isin(e+iDt ) ♦ P c ico3(e-wt) ♦ P a isln(e-wt) (74) 

4 4 -- 

w »ii • W^lcosCe+uthW^islnte+toO+Wjcicosfe-wtJ+W^isinCe-wt) (75) 

4 4 - - 

^aii ■ W a clc° s (0' , 'wt)*W J sisin(B*to>t)'*-W a clco3(0-wt)«-W a si3in(e-tot) (76) 
Substituting equations (73). (74), (75) and (76) Into equations (70), 
(71) and (72) and grouping like terms of sines and cosines (as shown In 
Appendix C) eliminates the time and theta dependency and yields twelve 
linear algebraic equations per cavity. The resulting system of 
equations for the i-th cavity can be stated: 

tA i-1 3 CXi- 1 ) ♦ [Ai] (Xi) ♦ [Ai*i3 (Xi*i) - a (Bi) ♦ b (Cl) (77) 

€ € 

where 

4 4 - - 4 4 - 

(x i-l) - (Psi-1. Pci-1. Psi-1, Pcl-1. W.si-i, W,ci-1 , W.si- 1 , 

w ici-1t W 2 si- 1 , W 2 si- 1 , W 2 si- 1 , W 2 ci-l) 

(Xi) - (Psi. Pci. Psi. Pci, W*si, W*d , W ,si , H~ci, W a si, W a ci , 
W asl. W a ci) 

4 4 - — 4 4 - 

(Xi+l) - (Psi + 1, Pci+1, Psi + 1, Pci-M. W lS i + i, W.ci + i, W,si + i , 

u - -*• + -- T 

"icl+1, W a si+ 1 , W 2 cl+1, Wjsi+i, Wjcl+l) 

The A matrices and column vectors B and C are given In Appendix C. To 

use equation (77) for the entire seal solution, a system matrix can be 

formed which is block tridiagonal in the A matrices. The size of this 

resultant matrix is (12NC X 12NC) since pressure and velocity 

perturbations at the inlet and the exit are assumed to be zero. This 

system is easily solved by various linear equation algorithns, and 

yields a solution of the form: 
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♦ a ♦ b ♦ 

■si ■ ~ Fasl * ~ Fbsi 

c c 

- a - b 

Pal “ - Fasi ♦ “ Fbsi 

t c ( 78 ) 

♦ a ♦ b ♦ 

Pci * " Faci * ~ Fbcl 

e c 

- a - b - 

p ci ■ ~ Facl ♦ “ Fbci 

e c 

Determination of Dynamic Coefficient 

The force-motion equations for a labyrinth seal are assumed to be 
of the form: 


- ck: .1 i;i •[:,:] i;i 


(79) 


The solution of equation (79) for the stiffness and damping 
coefficients is the objective of the current analysis. For the assumed 
elliptical orbit of equation (73). the X and Y components of 
displacement and velocity are defined as: 

X ■ a cosut X ■ -aii) sinut 

i 

Y ■ b si nut Y ” bu cosut 

Substituting these relations into equation (79) yields: 

F x ■ -Ka cosut - kb sinut ♦ Cau slnut - cbu cosut (80) 

Fy ■ -ka cosut - Kb sinut - cau slnut - Cbu cosut 

Redefining the forces, Fx and Fy, as: 

Fx * F xc cosut ♦ Fxs sinut (8 1 ) 

Fy “ FyQ cosut ♦ Fys sinut 

and substituting back into equation (78) yields the following 

relations: 


~F X c ■ Ka ♦ cbu “Fxs ■ ~Cau ♦ kb (82) 

“Fyc » ka ♦ Cbu ”Fys ■ Kb ♦ cau 

The X and Y components of force can be found by integrating the 


# 




« 


% 


pressure around the seal as follows: 
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NC 

2n 






Fx - -Rse 

l 

i-1 

/ 

0 

hi 

Li 

cose 

de 

(83) 


NC 

2* 






Fy - Rse 

l 

i-1 

J 

0 

hi 

Li 

sine 

de 

(84) 


Only one of these components needs to be expanded In order to determine 
the dynamic coefficients. For this analysis, the X component was 
chosen. Substituting equation (74) into ( 83 ) and Integrating yields: 

NC ♦ ♦ 

Fx * ~exR3 J Li [(Psi ~ Psi) sinwt ♦ (P C i 4 Pci) cosu»t] ( 85 ) 

1-1 

substituting from equations (78) and (80) into equation (85) and 

equating coefficients of sinwt and coswt yields: 

NC ♦ ♦ 

F xs ■ ~xRs J Li [a(Fasi ~ Fasi) 4 b (Fbsi ~ Fbsi)3 

i-1 (86) 

NC 

F xc - -xRs 1 Li Ca(Facl 4 Faci) 4 b (Fbci 4 Fbci)3 
1-1 

Equating the alternative definitions for F xa and Fxc provided by 
equations (82) and (86) and grouping like terms of the linearly 
Independent coefficients a and b yields the final solutions to the 
stiffness and damping coefficients: 



NC 

♦ - 

K 

- irR l 

(F a ci 4 Faci^ Lj 


1-1 



NC 

♦ - 

k 

- »R I 

( p bsi - Fbsi ^ Lj 


i-1 



-xRs NC 

4 — 

C 

• w l 

( F asl ~ Fasl^ Lj 


i-1 



xRs NC 

4 — 

c 

‘ “ i-f 

( F bci 4 Fbci ^ L i 


( 87 ) 



Data Requirements and Solution Procedure Summary 


The required input for the analysis presented is as follows: 

a) Reservoir pressure, temperature, and kinematic viscosity. 

b) Sump pressure. 

c) Gas constant and ratio of specific heats. 

d) Inlet circumferential velocity and rotor speed. 

e) Seal radius, radial clearance, tooth pitch, height and tip 
width. 

f) Rotor and stator friction coefficients (mr,nr ,ms,ns) . 

g) Number of teeth. 

In review, the solution procedure uses the following sequential steps: 

a) Determination of whether flow is choked or not using equations 
(63) and (66). 

b) The steady-state pressure distribution and leakage are found 
using equation (63)and/or (66). 

c) The steady-state circumferential velocity distribution is 
determined using equation (68). 

d) A system equation is formed for the first-order 
perturbation variables and solved using the cavity equation 
( 77 ). 

e) Results of this first-order perturbation solution, as defined 
in equations (78), are inserted into equation (87) to 
define the rotordynamlc coefficients. 


CHAPTER III 


TEST APPARATUS AND FACILITY 


TESTING APPROACH 


The testing method employed at the TAMU facility Is the same as 
that used by lino and Kaneko [HU]. An external hydraulic shaker is 
used to impart translatory motion to the rotating seal, while rotor 
motion relative to the stator and the reaction force components acting 
on the stator are measured. 

Figure 21 shows the manner in which the rotor could be positioned 
and oscillated in order to Identify the dynamic coefficients of the 
seal for small motion about an eccentric position, e 0 . Equation (1) is 
rewritten here 

C XXUo) K XY(e 0 n (Xi ^ r c xx(e 0 ) C XY(e 0 n (x) 

YX(e 0 ) K YY(co)J * y! Lc YX ( eo ) C YY ( eo J \y/ (88) 

First, harmonic horizontal motion of the rotor is assumed, where 

x “ e o ♦ A sin(fit) ♦ B cos(Ot) 

X - AO cos(Ot) - Bn sin (nt ) 

Y - 1 - 0 

This yields small motion parallel to the static eccentricity vector, 
where Q is the shaking frequency. In a similar fashion, the X and Y- 
dlrection force components can be expressed 



( 89 ) 


F X - FXS sin(nt) ♦ FXC cos(Qt) 
f Y - FYS Sln(nt) ♦ Fyc cos(Ot) 

Substituting these expressions into equation (88) and equating 
coefficients of constant, sine, and cosine terms yields the following 
four equations for the dynamic coefficients 
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F XS- K XX A - C XX B 
F XC “ K XX B 4 C XX A 

(90) 

F YS “ K XX A - Cy X B 
F YC * *YX B ♦ Cy X A. 

Solving this system of four equations in four unknowns defines the 
dynamic coefficients as 

K XX< E o> • < F XC B ♦ Fys A) / (A* ♦ B a ) 

Kyx(co> " < f YS * ♦ Fy C B) / (A a ♦ B a ) 

(91) 

c XX< e o) " ( F XC A - Fxs B) / Q(A a ♦ B a ) 

c YX^ E o) ■ (FYC A ~ Fys B) / Q(A a ♦ B a ) 

Therefore, by measuring the reaction forces due to known rotor 
motion, determining the Fourier coefficients (A,B,F X s»Fxc* F YS* F YC^ • and 
substituting into the above definitions, the indicated dynamic 
coefficients can be identified. If the rotor is shaken about a 
centered position (e o -0), the process is complete. Since the 
linearized model has skew-symmetric stiffness and damping matrices, all 
of the coefficients are identified. If, however, the rotor is shaken 
about an eccentric position as initially postulated, then it must be 
shaken vertically about that same point in order to complete the 
identification process. 

Assuming harmonic vertical motion of the rotor, as defined by 

X - e 0 , X - 0, 

Y • A sin(nt) ♦ B cos(Qt), and 
Y • Afi cos(ot) - Bn sln(Qt), 

yields oscillatory motion that Is perpendicular to the assumed static 


eccentricity vector. A similar process as before results in the 


coefficient definitions 


K YY(e 0 ) - (FXS A ♦ Fxc B) / (A a ♦ B a ) 

K XY(eo) - “(FYC B ♦ FYS A) / (A* ♦ B*) 

( 92 ) 

C YY(e 0 ) - (FXC A - FXS B) / Q(A* ♦ B a ) 
c XY(c 0 ) - (FYS B - Fyc A) / Q(A* ♦ B a ). 

All eight dynamic coefficients are thus determined by alternately 
shaking the rotor at one frequency ft In directions which are parallel 
and perpendicular to the static eccentricity vector. 

APPARATUS OVERVIEW 

Detailed design of the TAMU gas seal apparatus was carried out by 
J.B. Dressraan of the University of Louisville. It Is of the external 
shaker configuration, with the dynamlc-coeff lclent-identif lcation 

process described In the preceding section. 

Considering both the coefficient identification process and the 
analysis, some objectives for the design of the test apparatus are 
apparent. First, to determine the dynamic coefficients, the apparatus 
must provide for (a) the necessary rotor motion within the seal, and 
(b) measurement of the reaction-force components due to this 
motion. Secondly, it would be advantageous (for purposes of 
comparison) if the apparatus could provide the same variable seal 
parameters afforded by the analysis (i.e., pressures, seal geometry, 
rotor rotational speed, fluid prerotation, and rotor/stator surface 
roughness). With this capability, the influence of each Independent 
parameter could be examined and compared for correlation between 
theoretical predictions and experimental results. 

With these design objectives in mind, the discussion of the test 


apparatus Is presented In three sections. The first section, Test 
Hardware, describes how the various seal parameters are physically 
executed and controlled. For example, the manner In which the dynamic 
"shaking" motion of the seal rotor Is achieved and controlled is 
described in this section. The second section. Instrumentation, 
describes how these controlled parameters, such as rotor motion, are 
measured. Finally, the Data Acquisition and Reduction section explains 
how these measurements are used to provide the desired information. 

TEST HARDWARE 

This section deals only with the mechanical components and 
operation of the test apparatus. It provides answers for the following 
questions: 

1) How is the static position of the seal rotor controlled? 

2) How is the dynamic motion of the rotor executed and 

controlled? 

3) How is compressed air obtained and supplied to the apparatus, 
and how is the pressure ratio across the seal controlled? 

*0 How is the incoming air prerotated before it enters the seal? 

5) How are the seal rotor and stator mounted and replaced? 

6) How is the seal rotor driven (rotated)? 

Recalling the rotordynamic-coeff icient-identif lcatlon process 

described earlier, the external shaker method requires that the seal 
rotor be set in some static position and then oscillated about that 
point. The test apparatus meets those requirements by providing 
independent static and dynamic displacement control, which are 


described below. 


Static Displacement Control . 

The test apparatus Is designed to provide control over the static 
eccentricity position both horizontally and vertically within the seal. 
The rotor shaft is suspended pendulum-fashion from an upper, rigidly 
mounted pivot shaft, as shown in figures 22 and 23. This 
arrangement allows a slde-to-slde (horizontal) motion of the rotor, and 
a cam within the pivot shaft allows vertical positioning of the rotor. 

The cam which controls the vertical position of the rotor is 
driven by a remotely-operated DC gearhead motor, allowing accurate 
positioning of the rotor during testing. Horizontal positioning of the 
rotor is accomplished by a Zonlc hydraulic shaker head and master 
controller, which provide Independent static and dynamic displacement 
or force control. The shaker head is mounted on an I-beam support 
structure, and can supply up to 1**150 N (1000 lbf) static and kii50 N 
dynamic force at low frequencies. The. dynamic force decreases as 
frequency is increased. As illustrated in figure 22, the shaker head 
output shaft acts on the rotor shaft bearing housing, and works against 
a return spring mounted on the opposite side of the bearing housing. 
The return spring maintains contact between the shaker head shaft and 
the bearing housing, thereby preventing hammering or the shaker shaft 
and the resulting loss of control over the horizontal motion of the 
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Fig. 22 Components used for static and dynamic displacement of seal rotor. 
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Dynamic Displacement Control . 


The dynamic motion of the seal rotor within the stator is 
horizontal. In addition to controlling the static horizontal position 
of the rotor, the Zonlc shaker head moves the rotor through horizontal 
harmonic oscillations as the test is run. A Wavetek function generator 
provides the sinusoidal Input signal to the Zonlc controller, and both 
the amplitude and frequency of the rotor oscillations are controlled. 

Although the test-rig design provides for dynamic motion of the 
rotor only In the horizontal X-direction, all of the coefficients for 
either seal model (equation (1) or (2)) can still be determined. As 
figure 2 h shows, the required rotor motion perpendicular to the static 
eccentricity vector can be accomplished in an equivalent manner by 
statically displacing it the same amount (e 0 ) In the vertical direction 
and continuing to shake horizontally. 

In addition to providing control over the rotor's static position 
and dynamic motion, the test apparatus allows other seal parameters to 
be controlled independently, providing insight into the influence these 
parameters have on seal behavior. These parameters coincide with the 
variable input parameters for the analysis, and they include: 

1) pressure ratio across the seal, 

2) prerotation of the incoming fluid, 

3) seal configuration, and 

*0 rotor rotational speed. 



Pressure Ratio. 


The inlet air pressure and attendant mass flow rate through the 
seal are controlled by an electric-over-pneumatically actuated 
Hasonellan Camflex II flow control valve located upstream of the test 
section. An Ingersoll~Rand SSR-2000 single stage screw compressor 
rated at 3* mVmln 9 929 kPa (1200 scfm § 120 psig) provides compressed 
air, which Is then filtered and dried before entering a surge tank. 
Losses through the dryers, filters, and piping result In an actual 
maximum inlet pressure to the test section of approximately 825 kPa 
(105 psig) and a maximum flow rate of 10 mVmin (350 scfm). A 
foui inch inlet pipe from the surge tank supplies the test rig, and 
after passing through the seal, the air exhausts to atmosphere through 
a manifold with muffler. 

Inlet Circumferential Velocity Control . 

In order to determine the effects of fluid rotation on the 
rotordynamic coefficients, the test rig design also allows for 
prerotation of the incoming air as it enters the seal. This 
prerotation introduces a clrcunferential component to the air flow 
direction, and is accomplished by guide vanes which direct and 
accelerate the flow towards the annulus of the seal. Figure 
25 illustrates the vane configuration. Five sets of guide vanes are 
available; two rotate the flow in the direction of rotor rotation at 
different speeds , another introduces no fluid rotation, and two rotate 
the flow opposite the direction of rotor rotation at different speeds. 
The important difference between the vanes is the gap height, A. The 
vanes with a small gap height produce the highest inlet tangential 
velocity. 




Seal Configuration . 


The design of the test rig, figure 26, permits the Installation or 
various rotor/stator combinations. The stator Is supported in the test 
section housing by three Klstler quartz load cells In a trihedral 
configuration, as shown in figure 27. Different seal stator designs are 
obtained by the use of Inserts. The smooth and labyrinth inserts used 
for the .limn (.016 'in.) radial clearance seal tests are shown In figure 
28. The labyrinth rotor and the tooth detail are shown in figures 29 
and 30. Seals with different geometries (i.e., clearances, tapers, 
lengths) can be tested, as well as seals with different surface 
roughnesses. 

Rotational Speed . 

A Westlnghouse 50-hp variable-speed electric motor drives the 
rotor shaft through a belt-driven jackshaft arrangement. This shaft is 
supported by two sets of Torrington hollow-roller bearings [i»5] . These 
bearings are extremely precise, radially preloaded, and have a 
predictable and repeatable radial stiffness. The shaft bearings are 
lubricated by a positive-displacement gear-type oil pump. 

Different Jackshaft drive-pulleys can be fitted to provide up to a 
*1:1 speed increase from motor to rotor shaft, which would result in a 
rotor shaft speed range of 0-21 ,200 cpm. Previously, the maximum 
possible test speed was 8500 cpm. High bearing temperatures and the 
reduction or interference in the rotor-shaft fitment with increasing 
speed had served to limit shaft speed. These problems have been 
addressed by some specific design modifications which are discussed 


below. 
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Fig. 26 Cross-sectional view of test section 
showing rotor-shaft assembly. 
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Fig. 30 Detail of labyrinth tooth. 
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In the past, the seal rotor was press-fitted and secured axially 
by a bolt circle to the rotor shaft. As the running speed Is 
Increased, however, the inertia-induced diametral growth of the rotor 
exceeds the growth of the shaft. By increasing the interference in 
stationary rotor-shaft fit, a greater allowance for this growth 
difference has been provided. Figure 31 shows the present rotor-shaft 
design, a tapered rotor which is hydraulically expanded during 
installation. The rotor is inserted over the end of the tapered shaft 
and a large nut is used to pull the rotor onto the shaft. Fluid is 
pumped between the shaft and rotor, causing the rotor to expand. This 
separating force allows the rotor to be pulled onto the shaft until the 
desired interference fit is achieved. 

The problem of high bearing temperatures has been eliminated by 
replacing a roller-type thrust bearing and modifying the lubricant 
flow. A Torrlngton Hydraflex. thrust bearing, consisting of eight one- 
inch rubber-faced pads which are water lubricated, is now in place at 
the rear of the rotor. In addition, the lubricant for the Torrlngton 
hollow-roller bearings which support the shaft has been changed to 
light turbine oil with a maximum temperature of 270°F. The hollow- 
roller-bearing caps have been modified to direct the oil flow to the 
regions of heat buildup. These modifications are shown in figure 31. 

The final modification to allow operation of the TAMU gas seal 
test apparatus at high speeds was the installation of Koppers 
circumferential seals for the hollow-roller and thrust bearing 
lubrication systems. At 16,000 rpm, the surface speeds of the shaft 
and rotor (170 and 350 ft/sec, respectively) exceed the limits of lip 



Fig. 31 High speed rotor-shaft assembly. 
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seals, which had been used on the TAMU apparatus. The Koppers seals in 
figure 31 were designed for gas applications. The sealing mechanism is 
a segmented carbon seal ring. 

To conclude this discussion of the test hardware, two views of the 
complete test apparatus are included. Figure 32 shows the assembled 
rig, while an exploded view is provided in figure 33. 



Fig. 32 Test apparatus assembly. 










INSTRUMENTATION 


Having discussed the seal parameters that can be varied, and how 
the variations are Implemented, their measurement will now be 
described. The types of measurements which are made can be grouped into 
the following three categories: 

1) rotor motion, 

2) reaction-force measurements, and 

3) fluid flow measurements. 

These categories are described Individually in the sections that 
follow. 

Rotor Motion Measurements . 

The position of the seal rotor within the stator is monitored by 
four Bently-Nevada eddy-current proximity probes, mounted in the test 
section housing. These probes are located 90 degrees apart, and 
correspond to the X and Y- directions. The proximity probes are used 
to determine the static position and dynamic motion of the rotor, and 
their resolution is 0.0025 mm (0.1 mil). 

Reaction-Force Measurements . 

Reaction forces arise due to the motion of the seal rotor 
within the stator. The reaction forces (Fy, Fy) exerted on the stator 
are measured by the three Klstler quartz load cells which support the 
stator in the test section housing. When the rotor is shaken, vibration 
is transmitted to the test section housing, both through the thrust 
bearing and through the housing mounts. The acceleration of the 
housing and stator generates unwanted inertial H ma" forces which are 
sensed by the load cells, in addition to those pressure forces 


developed by the relative notion of the seal rotor and stator. For 
this reason, PCB piezoelectric accelerometers with integral amplifiers 
are mounted in the X and Y-dlrections on the stator, as shown in 
figure 27. These ^ccels allow a (stator mass) x (stator acceleration) 
subtraction to the forces (F x , Fy) indicated by the load cells. With 
this correction, which is described more fully in the next section, 
only the pressure forces due to relative seal motion are measured. 

Force measurement resolution is a function of the stator mass and 
the resolution of the load cells and accelerometers. Accelerometer 
resolution is 0.005 g, which must be multiplied by the stator mass in 
order to obtain an equivalent force resolution. The masses of the 
stators used in the test program reported here are 11.5 kg(25.3 lb) and 
11.0 kg (2 4. 2 lb), corresponding to the smooth and labyrinth stators, 
respectively. Hence, force resolution for the accelerometers is 0.560 
N (0.126 lb) and 0.538 N (0.121 lb), for each stator, respectively. 
Resolution or the load cells is 0.089 N(0.02 lb). Therefore, the 
resolution of the force measurement is limited by the accelerometers. 
With a stator with less mass, and/or accelerometers with greater 
sensitivity, force resolution could be improved. 

Fluid Flow Measurements . 

Fluid flow measurements Include the leakage (mass flow rate) of 
air through the seal, the pressure gradient along the seal axis, and 
the inlet fluid circumferential velocity. 

Leakage is measured with a Flow Measurement Systems Inc. turbine 
flowmeter located in the piping upstream of the test section. 
Resolution of the flowmeter is 0.0005 acf, and pressures and 
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temperatures up and downstream of the meter are measured for mass flow 
rate determination. 

For measurement of the axial pressure gradient, the stator has 
pressure taps drilled along the length of the seal in the axial 
direction. These pressures, as well as all others, are measured with a 
0-1. 03 1 * MPa (0-150 psig) Scanivalve differential-type pressure 
transducer through a *18 port, remotely-controlled Scanivalve model J 
scanner. Transducer resolution is 0.552 kPa (0.08 psl). Overall 
accuracy of the pressure measurements is limited by the resolution of 
the 12 bit A/D converter which can only resolve the pressure signal to 
^ 0.62 kPa (0.09 psl). Combined linearity and hysteresis error for the 
pressure transducer is 0.06J. 

In order to determine the clrcunferentlal velocity of the air as 
it enters the seal, the static pressure at the guide vane exit is 
measured. This pressure, in conjunction with the measured flowrate and 
inlet air temperature, is used to calculate a guide vane exit Mach 
number. A compressible flow continuity equation 

■ - p ex A ex Me X [(Y/RT t ) (1 ♦ (Y-1)M ex 2 / 2)] * /2 (93) 
is rearranged to provide a quadratic equation for Hex 

M ex 2 - {-1 ♦ 1 ♦ H((Y-1)/2Y) (m RT t /p ex A ex ) 2 ) * CY-1 ) (9*0 
where Y is the ratio of specific heats and R is the gas constant for 
air, Tt is the stagnation temperature of the air, p ex is the static 
pressure at the vane exit, and A ex is the total exit area of the guide 
vanes. Since all of the variables in the equation are either known or 
measured, the vane exit Mach number, and therefore the velocity, can be 
found. 
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In order to determine the circumferential component of this inlet 
velocity, a flow turning angle correction, in accordance with 
Cohen [A6], is employed. The correction has been developed from guide 
vane cascade tests, and accounts for the fact that the fluid generally 
is not turned through the full angle provided by the shape of the guide 
vanes. With this flow deviation angle calculation, the actual flow 
direction of the air leaving the vanes (and entering the seal) can be 
determined. Hence, the magnitude and direction of the inlet velocity 
is known, and the appropriate component is the measured inlet 
circumferential velocity. 

DATA ACQUISITION AND REDUCTION 

With the preceding explanations of how the seal parameters are 
varied, and how these parameters are measured, the discussion of how 
the raw data is processed and Implemented can begin. Data acquisition 
is directed from a Hewlett-Packard 9816 (16-blt) computer with disk 
drive and 9.8 megabyte hard disk. The computer controls an H-P 69MOB 
multi programmer which has 12-bit A/D and D/A converter boards and 
transfers control commands to and test data from the instrumentation. 

As was previously stated, the major data groups are seal 
motion/reaction force data and fluid flow data. The motion/reaction 
force data are used for dynamic coefficient identification. The 
hardware Involved Includes the load cells, accelerometers, X-direction 
motion probe, a Sensotec analog filter unit, a tuneable bandpass 
filter, and the A/D converter. The operation of these components is 
illustrated in figure 3*1, and their outputs are used in a serial 
sampling scheme which provides the computer with the desired data for 



Pig. 34 Signal conditioning schematic for data acquisition* 


















reduction. Recalling the discussion of the reaction force measurements 
In the preceding section, a (stator mass) x (stator acceleration) 
subtraction from the Indicated load cell forces is necessitated due to 
vibration of the stator and test section housing. This subtraction is 
performed with an analog circuit, and results in corrected Fx and Fy 
force components due to relative seal motion. The forced oscillatory 
shaking motion of the seal rotor is the key to the operation of the 
serial synchronous sampling (SSS) routine which is employed. The 
frequency of the rotor oscillation is set by a function generator, and 
rotor motion is sensed by the X-direction motion probe. The motion 
signal is filtered by the narrow bandpass filter, and is used as a 
trigger signal for the SSS routine. Upon the operator’s command, the 
SSS routine is enabled, and the next positive-to- negative crossing of 
the filtered motion signal triggers a quartz crystal clock/timer. Ten 
cycles of the corrected Fy(t) signal are sampled, at a rate of 100 
samples/cycle. The second positive-to-negative crossing of the 
filtered motion signal triggers the timer and initiates the sampling of 
ten cycles of the Fy(t) signal. Finally, the third positlve-to-negatlve 
crossing triggers the timer again, and ten cycles of the corrected X(t) 
signal are sampled. Thus, at every test condition, 1000 data points 
are obtained for FyCtj ),Fy(t j ), and X(ty), and the data arrays are 
stored in computer memory. 

Some important points need to be stressed concerning this 
force/motion data acquisition. First, the bandpass filter is used only 
to provide a steady signal to trigger the timer/clock. Any modulation 
of the motion signal due to rotor runout is eliminated by this filter, 


as long as the rotational frequency and shaking frequency are 
adequately separated, and the shaking frequencies are selected to 
provide adequate separation with running speeds. However, the rotor 
notion and corrected force signals which are sampled and captured for 
coefficient identification are filtered only by a low-pass filter (500 
Hz cutoff), and the effects of runout as well as shaking motion are 
present In the recorded data. A second point worth noting Is that the 
sample rate Is directly dependent on the shaking frequency. As the 
shaking frequency Is increased, the sample rate (samples /second) also 
increases. In order to get the desired 100 samples /cycle, shaking 
frequencies must be chosen to correspond to discrete sample rates which 
are available. Hence, the frequency at which the rotor is shaken Is 
carefully chosen to provide the desired sampling rate and a steady 
trigger signal. The uncertainty in the shaking frequency is 0.13 Hz 
for the 74.6 Hz case. 

Most of the fluid flow data are used for the input parameters 
required by the analysis. The upstream (reservoir) pressure and 
temperature, downstream (sump) pressure, and the inlet circumferential 
velocity (determined as outlined earlier) are provided directly. The 
friction-factor values of the rotor and stator are supplied in the form 
of coefficients, which are obtained from the pressure distribution data 
for the smooth annular seals, see Nicks [473 and Nelson et al. [48], 
and are assumed to be the same for the labyrinth surfaces. 

PROCEDURE 

At the start of each day’s testing, the force, pressure, and 
flowmeter systems are calibrated. The total system, from transducer to 


computer, is calibrated for each of these variables. The force system 
calibration utilizes a system of pulleys and known weights applied in 
the X and Y-dlrections. An air-operated dead-weight pressure tester is 
used for pressure system calibration, and flowmeter system calibration 
is achieved with an internal precision clock which simulates a known 
flowrate. 

All of the tests performed to date have been made with the rotor 
executing small motion about a centered position. A typical test 
begins by centering the seal rotor in the stator with the static 
capability of the Zonlc hydraulic shaker, starting the airflow through 
the seal, setting the rotational speed of the rotor, and then beginning 
the shaking motion of the rotor. Data points are taken at rotational 
speeds of 3000 , 6000 , 9500, 13000 and 16000 cpm with a tolerance of 
10 cpm. At each rotational speed, data points are taken at pressures 
of 3.08 bar (30 pslg), H.H6 bar (50 psig), 5.8 1 * bar (70 pslg), 7.22 bar 
(90 psig), and 8.25 bar (105 pslg), as measured upstream of the 
flowmeter with a tolerance of ± 0.069 bar (1.0 pslg). For each test 
case (i.e., one particular running speed, shaking frequency, inlet 
pressure, and prerotation condition), the measured leakage, 
rotordynamlc coefficients, and axial pressure distribution are 
determined and recorded. 

This test sequence is followed for each of two different shaking 
frequencies, and for five inlet swirl directions. Therefore, twenty- 
five data points are taken per test with a total of ten tests per seal 
for shaking about the centered position. Shaking a seal about an 
eccentric position would require more tests. 


CHAPTER IV 


TEST RESULTS: INTRODUCTION 

The results reported here are from tests of six "see-through" 
labyrinth seals, three with teeth on the rotor and three with teeth on 
the stator, each with different radial clearances. Tables 3, A and 5 
show the pertlnant data for each seal configuration. For the remainder 
of this report, the seals will be referred to as seal 1, seal 2, and 
seal 3, as given in table 5, in addition to their respective 
configuration. 

The test program had the following objectives: 

1) Acquire leakage, stiffness, and damping coefficients as a 
function of rotor speed, pressure drop, and inlet circumferential 
velocity for three teeth-on-rotor and three teeth-on-stator 
labyrinth seals with different radial clearances. 

2) Compare the effect of varying the radial seal clearance 
on the experimentally determined rotordynamic coefficients. 

3) Compare test results to the predictions of the new analysis 
presented in this report. 

When shaking about the centered position, the test apparatus can 
be used to control the rotor speed, reservoir pressure (i.e. supply 
pressure), circumferential velocity of the inlet air, and the 
frequency and amplitude of translatory rotor motion. Two shake 
frequencies, 56.8 and 74.6 Hz, were used during testing with 
essentially the same results. The results plotted here were obtained 


Table 3. Test stator specifications. 


Seal 1 

Diameter : 
upstream 
downstream 

Material: 


Smooth Stator 


Labyrinth Stator 


15.197 cm (5.983 in) 
15.197 cm (5.983 in) 

aluminum 


15.202 cm (5.985 in) 
15.202 cm (5.985 in) 

brass 


Seal 2 

Diameter: 

upstream 

downstream 

Material : 


15.217 cm (5.991 in) 

15.217 cm (5.991 in) 

aluminum 


15.217 cm (5.991 in) 

15.217 cm (5.991 in) 

brass 


Seal 3 

Diameter : 
upstream 
downstream 

Material: 


15.245 cm (6.002 in) 
15.247 cm (6.003 in) 

brass 


15.237 cm (5.999 in) 
15.237 cm (5.999 in) 

brass 


Table 4. Test rotor specifications. 


Seal 

Diameter: 

upstream 

downstream 

Material: 


Labyrinth Rotor 


Smooth Rotor 


1.2,3 


1.2,3 


15.136 cm (5.959 in) 

15.136 cm (5.959 in) 

304 stainless steel 


15.136 cm (5.959 in) 

15.136 cm (5.959 in) 

304 stainless steel 


Table 5. Test seal specifications 


Seal 1 

Radial Clearance: 
upstream 
downstream 

Seal Length: 

Number of teeth: 

Seal 2 

Radial Clearance: 
upstream 
downstream 

Seal Length: 

Number of teeth: 

Seal 3 

Radial Clearance: 
upstream 
downstream 

Seal Length: 

Number of teeth: 


Teeth-On-Rotor 

0.3048 cm (0.012 in) 
0.3048 cm (0.012 in) 

5.080 cm (2.000 in) 

16 

0.4064 cm (0.016 In) 
0.4064 cm (0.016 in) 

5.080 cm (2.000 in) 

16 


Teeth-On-Stator 

0.3302 cm (0.013 In) 
0.3302 cm (0.013 in) 

5.080 cm (2.000 In) 

16 

0.4064 cm (0.016 in) 
0.4064 cm (0.016 in) 

5.080 cm (2.000 In) 

16 

0.5080 cm (0.020 In) 
0.5080 cm (0.020 in) 

5.080 cm (2.000 In) 

16 


0.5461 cm (0.0215 in) 
0.5588 cm (0.022 in) 

5.080 cm (2.000 in) 

16 


Table 6. Definition of symbols used in figures. 


Pressures 

1 3»08 bar 

2 - 4.46 bar 

3 - 5.84 bar 

4 - 7.22 bar 
5-8.25 bar 


Rotor speeds 

1 - 3000 cpm 

2 - 6000 cpm 

3 “ 9500 cpm 

4 - 13000 cpm 

5 - 16000 cpm 


Inlet circumferential velocities 

1 <- High velocity against rotation 

2 - Low velocity against rotation 

3 “ Zero circumferential velocity 

4 - Low velocity with rotation 

5 - High velocity with rotation 


The pressure for each test is set at the flowmeter of Figure 12. 


by shaking at 7*1.6 Hz at an amplitude between 3 and *1 mils. The actual 
test points for each of the other three independent variables are shown 
in Table 6. 

Figures 35-37 show the inlet circumferential velocity of the air 

(U e ) for the configurations described in table 6 for the seals reported 
on here. The equation for U e is 

u 6 ■ m sin A / p A v 

where m is the fluid mass flow rate, p is the fluid density, A v i s the 
exit area of the fluid turning vanes, and A is the fluid swirl angle at 
the turning vanes exit as measured from the axial direction. The 
method used to determine A is described in the TEST APPARATUS AND 

I 

FACILITY chapter of* this report. Negative circumferential velocities 
represent velocities opposed to the direction of rotor rotation. 
Positive velocities are in the direction of rotor rotation. Note that 
curve 3 (representing zero inlet circumferential velocity) lies on the 
horizontal axis in each figure. The inlet circumferential velocity 
ratio, the ratio of inlet circumferential velocity to rotor surface 
velocity, ranged from about -6 to about 6. When reviewing the 

following figures, table 6 and figures 35~37 should be consulted for 
the definitions of symbols used. 

NORMALIZED PARAMETERS 

Before the tests described herein were performed, the TAMU gas 
seal test apparatus was modified as described in the TEST APPARATUS AND 
FACILITY chapter to allow operation at running speeds up to 16,000 cpm. 
As expected, subsequent tests revealed a dependence of the rotor 
diameter on running speed due to inertia and thermal effects. The 









Table 7. Oowth of rotor with rotational speed 


Rotor speed 
(rpm) 


Diametrical growth 
(are) (Inches x 1000) 


3000 

6000 

9500 

13000 

16000 


0.01 

0.11 

0.02 

0.7 

0.03 

1.2 

0.05 

1.9 

0.11 



Table 8. Normalized coefficients. 


H 0 

K - K (nondlm) 

DL(AP) 

H 0 

C - C ° (sec) 
DL(AP) 


- stiffness (N/mm) 

- seal exit clearance (mm) 

- seal length (m) 


C » damping (N sec/mm) 

D ■ seal diameter (m) 

AP - pressure drop across 
seal (N/m 1 ) 


rotor growth data, shown in table 7» were obtained from eddy 
current motion probes positioned at the mldspan of the seal. Thus, as 
the rotor turns faster, the forces in the seal are affected not only by 
the Increased surface speed of the rotor (drag) but also by a change In 
clearance (friction factor). See table 8 for the definitions of the 
normalized parameters. Theoretically, normalization would collapse the 
data and make the presentation simpler and more straight forward. 
However, this is not the case with the labyrinth seals tested in this 
study. Figure 38 shows a comparison of dimensional and nondimensional 
direct stiffness versus rotor speed for the inlet pressure set of table 
6. The data did not collapse to a single curve and shows increased 
irregularity . 

Figure 39 shows a comparison of normalized and dimensional direct 
damping for a teeth-on-stator labyrinth seal versus clearance for the 
inlet circumferential velocity set of table 6. The normalized results 
lead one to believe that the direct damping coefficient increases as 
clearance increases. However, the dimensional results show that the 
direct damping coefficient decreases as clearances Increases. To avoid 
this type of confusion, the rotordynamic coefficients presented in this 
study have not been normalized. 

RELATIVE UNCERTAINTY 

Before the test results are given, a statement about the 
experimental uncertainty is needed. The method used is that described 
by Holman [N93 for estimating the uncertainty in a calculated result 
based on the uncertainties in primary measurements. The uncertainty wr 
in a result R which is a function of n primary measurements 
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Pig. 38 A comparison of dimensional and nondlmenslonal direct 
stiffness coefficients. See table 8 for definitions of 
normalized parameters. 
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Fig. 39 A comparison of dimensional and normalized direct 
damping coefficients. See table 8 for definitions of 
normalized parameters. 
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In this case, the rotordynamic coefficients are calculated using 
equation (91). The primary measurements are forces, displacements, and 
frequency. The uncertainty In these measurements on the TAMU test 
apparatus are 0.89 N (0.2 lb), 0.0013 an (0.05 mils), and 0.13 Hz, 
respectively. For the six seals tested, the estimated uncertainty In 
the stiffness and damping coefficients were 7 N/mm (lio lb/ln) and 
0.0875 N-s/mm (0.5 lb-s/in), respectively. The uncertainty In the 
cross-coupled damping coefficients were of the same order of magnitude 
as the coefficients themselves. Since the uncertainties in the cross- 
coupled-damping values were so high, and since the cross-coupled- 
damplng forces are of minor significance compared to the other damping 
and stiffness forces, comparisons of the cross-coupled-damping 
coefficients have been omitted from this report. 

SELECTION OF REPORT DATA 

For each of the six seals tested, there were 125 test points for 
leakage, direct and cross-coupled stiffness, and direct damping at the 
7*1.6 Hz shake frequency. Generally, a ranking of the three independent 
variables of the test apparatus in order of the relative effect on the 
rotordynamic coefficients of a seal is: inlet circumferential velocity, 
pressure ratio, running speed. The previous report of Scharrer [22] 
thoroughly catalogued the results for the effects of pressure ratio, 
rotor speed up to 8000 cpm and inlet circumferential velocity on the 
rotordynamic coefficients. Since the rotor speed capability of the 



test apparatus has been changed, results which show new information 
concerning rotor speed will be reviewed. Figures in the next two 
chapters show the dependence of leakage and rotordynamic coefficients 
on radial seal clearance for Inlet swirl conditions of table 6. 
Figures in Appendix D show additional information on leakage, pressure 
gradients and rotordynamic coefficients. Generally, solid lines in a 
figure represent experimental results, and broken lines represent the 
predictions of the new analysis presented in this report. These figures 
will be used to compare the effect of radial seal clearance on seal 
performance, and to evaluate the new analysis presented in this report. 


CHAPTER V 


TEST RESULTS: RELATIVE PERFORMANCE OF SEALS 

This evaluation of the effect on seal performance of varying the 
radial seal clearance requires frequent use of the information in 
table 6 and figures 35-37. It might seem obvious, since this report 
evaluates the effect of radial seal clearance, that the data should 
be presented as a function of clearance (clearance being on the x- 
axls). However, since the inlet circumferential velocity is directly 
dependent on the seal leakage and the seals leak at different rates 
due to differing cross-sectional areas, inlet circumferential 
velocity 5 (swirl 5) for seal 1 is less than those for seals 2 and 3. 
This is a problem because the rotordynamic coefficients are very 
sensitive to the inlet circumferential velocity. Therefore, the 
dynamic data will be presented as a function of inlet circumferential 
velocity at one inlet pressure and one rotor speed. Comparisons of 
the leakage, direct stiffness, cross-coupled stiffness, direct 
damping, and stability of the six seals follow. 

LEAKAGE 

The flow rate of air through each seal was measured with a 

turbine flowmeter located in the piping upstream of the test section 

« 

(see figure 32). Figures MCM3 show seal leakage as a function of 
radial seal clearance for the inlet circumferential velocity set of 
table 6 for teeth-on-rotor and teeth-on-stator labyrinth seals, 
respectively. The plot on the left side of the page is for the 
teeth-on-rotor seal and the one for the teeth-on-stator seal is on 
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Pig. *0 Leakage versus radial seal clearance at an Inlet pressure of 
3.08 bar and rotor speed of 3000 cpm. See table 6 for Inlet 
swirl definitions. Teeth-on-rotor (left), teeth^on-stator 
(right) . 
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rig. m Leakage veraua radial seal clearance at an Inlet pressure of 
3.08 bar and rotor speed of 16000 cpn. See table 6 for inlet 

swirl definitions. Teeth-on~rotor (left) , teeth-on-stator 
(rignt ) . 
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Fig. 42 Leakage versus radial seal clearance at an Inlet pressure of 
8.25 bar and rotor speed of 3000 epa. See table 6 for Inlet 
swirl definitions. Teeth-on-rotor (left), teeth-on-stator 
(right) . 
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the right. This convention will be followed for the remainder of the 
presentation. A comparison of the leakage of the six seals reveals, as 
expected, greater leakage occurs for greater radial clearances. The 
difference In the shapes of the curves, between the two seal types, can 
be attributed to the difference in the performance characteristics of 
the two seals. 

DIRECT STIFFNESS 

Figures M-47 show the dimensional direct stiffness versus Inlet 

circumferential velocity ratio for the clearances of table 5. These 

plots show that the direct stiffness decreases in magnitude as 

clearance Increases. One would expect zero direct stiffness values at 

sufficiently large clearances. Figures 48-50 show the dimensional 

direct stiffness versus rotor speed for the pressure ratios of table 6 

and inlet circumferential velocity 3- The figures show that direct 

stiffness becomes increasingly negative as rotor speed Increases, for 

the teeth-on-rotor seal, and is unchanged for the teeth-on-stator seal. 

This effect was not noticeable in the results from the low speed test 

/ 

rig and could be a result of clearance change due to rotor growth. The 
dimensionless direct stiffness coefficient, defined in table 8, removes 
the effect of clearance change due to rotor growth from the plot. 
Figure 51 shows the dimensionless direct stiffness versus rotor speed 
for seal 1 (minimum clearance seal) at the pressure ratios of table 6 
and inlet circumferential velocity 3. The figure shows that the dimen- 
sionless direct stiffness increases in magnitude as rotor speed 
increases, for the teeth-on-rotor seal, and is inconclusive for the 
teeth-on-stator seal. Associated direct stiffness plots can be found in 
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Fig. ** Direct stiffness versus inlet circumferential velocity ratio 
at an inlet pressure of 3.08 bar and rotor speed of 3000 cpm. 

See table 5 for seal clearance definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Direct stiffness versus inlet circumferential velocity ratio 

* P p ® s * ur « of 3.08 bar and rotor speed of 16000 
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rotor (left), teeth-on-stator (right). 
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Fig. *6 Direct stiffness versus Inlet circumferential velocity ratio 
at an Inlet pressure of 8.25 bar and rotor speed of 3000 cpm. 
See tabie 5 for seal clearance definitions. Teetb-on-rotor 
(left), teeth-on-stator (right). 
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Fig. *7 Direct stiffness versus Inlet circumferential velocity ratio 
at an inlet pressure of 8.25 bar and rotor speed of 16000 
. C P®* See table 5 for seal clearance definitions. Teeth-on- 
rotor (left), teeth-on-stator (right). 
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Fig. *8 Direct stiffness versus rotor speed for seal 1 and inlet 
circumferential velocity 3» See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. $9 Direct stiffness versus rotor speed for seal 2 and Inlet 
circumferential velocity 3- See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Direct stiffness versus rotor speed for seal 3 and Inlet 

2i^!^f? rentla i Veloclty 3. See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. 51 Dimensionless direct stiffness versus rotor speed for seal 1 
and Inlet circumferential velocity 3. See table 6 for 
pressure ratio definitions. Teeth-on-rotor (left), teeth-on- 
stator (right). 
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CROSS-COUPLED STIFFNESS 


Figures 52-55 show cross-coupled stiffness versus inlet 
circumferential velocity ratio for the radial clearances of table 5. 
The plots show, in all cases, that the cross-coupled stiffness is a 
somewhat linear function of inlet circumferential velocity ratio. The 
figures do not show any consistent trend with respect to radial seal 
clearance, for either teeth-on-rotor or teeth-on-stator seals. Figures 
56-58 show cross-coupled stiffness versus rotor speed for the pressure 
ratios of table 6. The figures show that cross-coupled stiffness 
increases with increasing rotor speed, for the teeth-on-rotor seal, and 
decreases with increasing rotor speed for the teeth-on-stator seal. 
This effect was not evident in the results from the low speed test rig. 
The dimensionless cross-coupled stiffness coefficient, defined in table 
8, removes the effect of change of clearance due to rotor growth. 
Figure 59 shows dimensionless cross-coupled stiffness versus rotor 
speed for seal 1 for the pressure ratios of table 6 and inlet 
circumferential velocity 5. The figure shows that cross-coupled 
stiffness increases for increasing rotor speed, for the teeth-on-rotor 
seal, and decreases with Increasing rotor speed for the teeth-on-stator 
seal. The decrease in cross-coupled stiffness with rotor speed, for a 
teeth-on-stator seal, was also evident in tests of an 11 cavity seal 
for Sulzer [50] and in the steam tests of 1-3 cavity seals by Hisa et 
al. [83. Associated cross-coupled stiffness plots can be found in 
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rig. 52 Ooaa-coupled atiffneaa veraua Inlet circumferential velocity 
ratio at an Inlet preaaure of 3-08 bar and rotor apeed of 
3000 cpa. See table 5 for aeal clearance deflnltlona. 
Teeth-on-rotor (left), teeth-on-atator (right). 
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rig. 53 Ooaa-coupled atiffneaa veraua inlet circumferential velocity 

at an lnlet P resaur « of 3.08 bar and rotor speed of 
16000 cpm. See table 5 for seal clearance definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. 5* Cross-coupled stiffness versus Inlet circumferential velocity 
ratio at an Inlet pressure of 8.25 bar and rotor speed of 
3000 cpm. See table 5 for seal clearance definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. 55 Cross-coupled stiffness versus Inlet circumferential velocity 
ratio at an Inlet pressure of 8.25 bar and rotor speed of 
16000 cpm. See table 5 for seal clearance definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 


113 


Kxy CN/mm) 


SERL 1 SEAL 1 

SWIRL 3 SWIRL 3 



2000 G000 10000 14 000 2000 G000 10000 14000 


4000 8000 12000 1G000 4000 8000 12000 16000 

ROTOR SPEED (CPM) ROTOR SPEED (CPM) 

Fig. 56 Croaa-coupled atiffneaa veraua rotor apeed for aeal 1 and 
Inlet circumferential velocity 3. See table 6 for preaaure 
ratio deflnltlona. Teeth-on-rotor (left), teeth-on-etator 
(right). 
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Fig. 57 Cross~coupled stiffness versus rotor speed for seal 2 and 
inlet circumferential velocity 3. See table 6 for pressure 
ratio definitions. Teeth-on-rotor (left), teeth-on-stator 
(right) . 
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Fig. 58 Croas~coupled stiffness versus rotor speed for seal 3 and 
inlet circumferential velocity 3. See table 6 for pressure 
ratio definitions. Teeth-on-rotor (left), teeth-on-stator 
(right). 
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Fig. 59 Dimensionless cross-coupled stiffness versus rotor speed for 
seal l and Inlet circumferential velocity 5. See table 6 for 
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DIRECT DAMPING 


Figures 60-63 show the direct damping versus Inlet clrcumrerentlal 
velocity ratio for the radial seal clearances of table *1 . The data 
show that the direct damping for a teeth-on-rotor seal increases as 
clearance is increased while the direct damping for a teeth-on-stator 
seal decreases as clearance is increased. Figures 6*1-66 show direct 
damping versus rotor speed for the pressure ratios of table 6 and inlet 
circumferential velocity 3. These figures show that the direct damping 
increases slightly as rotor speed increases, for teeth-on-rotor seals, 
and decreases slightly as rotor speed increases for teeth-on-stator 
seals. These results are deceiving, since direct damping is very sen- 
sitive to clearance change. The normalized direct damping coefficient, 

defined in table 8, removes the effects of clearance change due to 

\ 

rotor growth. Figure 67 shows normalized direct damping versus rotor 
speed for seal 1 for the pressure ratios of table 6 and inlet 
circumferential velocity 3. This figure shows that direct damping 
decreases as rotor speed increases for both teeth-on-rotor and teeth- 
on-stator seals. The result for normalized direct damping versus rotor 
speed for the teeth-on-rotor seal is inconsistent with the dimensional 
data. If dimensional direct damping, for a teeth-on-rotor seal, 
decreases as clearance decreases and increases as rotor speed increases 
then the normalized value should increase as rotor speed increases 
because the seal clearance decreases as rotor speed increases. This 
inconsistency is not readily explainable. Associated direct damping 
plots can be found in Appendix D. 
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Fig. 60 Direct damping versus inlet circumferential velocity ratio 

1 ”J et P p e3aure of 3.08 bar and rotor speed of 3000 cpm. 
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Pig. 61 Direct damping versus Inlet circumferential velocity ratio 
at an Inlet pressure of 3.08 bar and rotor speed of 16000 cpm. 

See table 5 for seal clearance definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Fig. 63 Direct damping versus Inlet circumferential velocity ratio 
at an Inlet pressure of 8.25 bar and rotor speed of 16000 cpm. 

/?*~? able 5 for 8eal clearance definitions. Teeth-on^rotor 
(left), teeth-on-stator (right). 
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Fig. 65 Direct damping versus rotor speed for seal 2 and Inlet 
circumferential velocity 3. See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. 66 Direct damping versus rotor 

circumferential velocity 3. See 
definitions. Teeth-on-rotor (left 
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Fig. 67 Normalized direct damping versus rotor speed for seal 1 and 
Inlet circumferential velocity 3. See table 6 for pressure 
ratio definitions. Teeth-on-rotor (left), teeth-on-stator 
(right) . 
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STABILITY ANALYSIS 

One further parameter of comparison among the test seals is the 
dimensionless whirl frequency ratio. To understand the value of this 
parameter , consider a rotor In a circular orbit of amplitude A and 
frequency u (Fig. 68). The X and Y components of force In the seal 

model of equation (79) may be resolved into radial and tangential 
forces 

F r ■ Fx cos cut ♦ Fy sin ut 

Ft ■ “Fx sin wt ♦ Fy cos ut 

Expressing the rotor motion as 

» 

X ■ A cos wt X « -Aw sin cut 

If * A sin wt Y « Aw cos cot 

and using equation (79), the resultant radial and tangential forces are 
Illustrated in the figure and are defined by 

“Fp/A m K ♦ cu 

F t /A “ k - Coi 

If F t /A is a positive quantity, the tangential force is destabilizing 
since It supports the whirling motion of a forward whirling rotor. 

Conversely, of F t /A is negative, it opposes the whirling motion of a 
forward whirling rotor, and is therefore stabilizing. The whirl 
frequency ratio is defined by 

Whirl frequency ratio » k/C« . 

From the above discussion, if the whirl ratio is less than one, the 
tangential force on the rotor is stabilizing. A minimum value of the 
whirl frequency ratio is optimum for stability. 








Figures 69 and 70 show the whirl frequency ratios at a running 
speed of 16000 cpm and a shake frequency of 7*».6 Hz. For teeth-on-rotor 
seals, the figures show that as clearance Increases the seal becomes 
more stable. For teeth-on-stator seals the opposite is true; as 
clearance Increases the seal becomes less stable, for the positive 
inlet circumferential velocity case. The figures also show that the 
teeth-on-stator seals are more stable than the teeth-on-rotor seals for 
positive inlet circumferential velocity ratio, as was found previously 
[63. Figure 71 shows the whirl frequency ratio versus rotor speed for 
the seals of table 5. The figure shows that as speed increases, both 
the teeth-on-rotor and teeth-on-stator seals become more stable. 
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Fig. 69 Whirl frequency ratio versus Inlet circunferentlal velocity 
at inlet Pressure of 3.08 bar and rotor speed of 

16000 cpra. Shake frequency 7*1.6 Hz. See table 5 for seal 
clearance definitions. Teeth-on-rotor (left), teeth-on-stator 
riant ) . 
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Fig. 70 Whirl frequency ratio versus Inlet clrcuaferential velocity 
ratio at an Inlet pressure of 8.25 bar and rotor speed of 
16000 cpn. Shake frequency 7*1.6 H*. See table 5 for seal 
clearance definitions. Teeth-on-rotor (left), teeth-on-stator 
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Fig. 71 Whirl frequency ratio versus rotor speed for the seals of 
table 5 at an inlet pressure of 8.25 bar and inlet 
circimferential velocity 5. 


CHAPTER VI 


TEST RESULTS: COMPARISON TO THEORETICAL PREDICTIONS 

In this chapter , the experimental results from the tests of three 
teeth-onrrotor and three teeth-on-stator labyrinth seals are compared 
to the new analysis presented In the Geometric Boundary Approach 
section in the THEORETICAL DEVELOPMENT chapter of this report. The 
seals tested are described in tables 3~5. Tables 5 and 6 and figures 
35-37 define the symbols used in the figures. Generally, the solid 
lines are the experimental points and the broken lines are the 
predictions . 

STATIC RESULTS 

Before proceeding with the comparison to the theory, some 
necessary input parameters to the model must be given. Table 9 shows 
the variables used as input to the program for the comparisons shown 
here. The temperature given was fairly constant for all of the tests. 
The viscosity was calculated for each case using Sutherland's formula 
[51], The pressure gradients for the five rotor speeds of table 6 are 
shown in figure 72 for a single inlet pressure and inlet 
circumferential velocity. The curves show that the pressure gradient 
has little or no sensitivity to rotor speed. Any slight differences In 
the curves are due to variations in the actual points taken. 
Therefore, only one rotor speed will be used for comparison of the 
pressure gradients. 

Figures 73~75 show a comparison of the experimental and 
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Table 9. Input parameters for seal program. 

Inlet temperature 300K 

Ratio of sp. heats l.ii 

Gas constant (air) 287.06 J/(kgK) 

Compressibility factor 1.0 

Rotor friction exp.(mr) - 0.25 

Rotor friction const. (nr) 0.079 

Stator friction exp. (ms) -0.25 

Stator friction const. (ns) 0.079 

Number of teeth 16 
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Fig. 73 A comparison of experimental and theoretical pressure 
gradients of seal 1 for a rotor speed of 3000 cpm and inlet 
circumferential velocity 3. See table 6 for inlet pressure 
definitions. Teeth-on^otor (left), teeth-on-stator (right). 
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Fig. 75 A comparison of experimental and theoretical pressure 
gradients of seal 3 for a rotor speed of 3000 cpra and Inlet 
circumferential velocity 3- See table 6 for inlet pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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theoretical pressure gradients for the no-lnlet-clrcumferential- 
velocity case. The figures show that the theory underpredicts the 
cavity pressures for the teeth-on-rotor case and overpredicts the 
cavity pressures for the teeth-on-stator case. This difference is due 

to the difference in the inlet losses for the two seal types. The 
theory accurately predicts the inlet loss for the teeth-on-rotor case. 
However, the teeth-on-stator seal has a much larger inlet loss followed 
by a pressure recovery. This positive pressure recovery cannot be 
modelled by a simple leakage equation. The remainder of the pressure 
gradient comparison plots can be found in Appendix D. 

Figures 76-78 show a comparison of experimental and theoretical 
leakage versus inlet circumferential velocity ratio for the inlet 
pressure set of table 6. The plots show that the theory underpredicts 
the leakage for both seal types by about 25%. This is much worse than 
the 5% error for ,the theory of Childs and Scharrer [18]. This 
difference is due to the change in the equation for the kinetic energy 
carryover coefficient, y a . The change in the coefficient was made in 
order to obtain a local equation which would yield a clearance 
perturbation. The former coefficient was a global equation and could 
not be perturbed. The contribution of this coefficient to the first- 
order equations and the subsequent solution is very substantial. In 
effect, the leakage calculation was sacrificed In order to improve the 
calculation of the dynamic coefficients. 

DYNAMIC RESULTS 

The experimental and theoretical results to be compared include 
the direct and cross-coupled stiffness and direct damping coefficients. 
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Fig. 76 A comparison of experimental and theoretical leakage 
versus inlet circumferential velocity ratio for seal 1 at a 
rotor speed of 3000 cpm. See table 6 for Inlet pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fl8 ‘ 78 ,° r experimental and theoretical leakage 

versus inlet circumferential velocity ratio for seal 3 at a 

def?nff?« eed r f J 000 C?m ' See uble 6 for lnl ®t pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right) 



A cross-coupled damping comparison has been omitted because or the 
uncertainty present in the experimental values (see the Relative 
Uncertainty section). Of the remaining three coefficients, the direct 
stiffness comparison will be presented first, and the direct damping 
comparison will be last. 

Direct stiffness 

Figures 7SN82 show a comparison of experimental and theoretical 
direct stiffness versus inlet circumferential velocity ratio for the 
seals defined in table 5. The figures show that the theory correctly 
predicts a decrease in direct stiffness as clearance increases, for 
both seal types. The figures also show that for both inlet pressures, 
the theory overpredicts the direct stiffness at low rotor speeds and 
underpredicts the direct stiffness at high rotor speeds. This trend Is 
made clearer by figures 83 - 85 . Figures 83-85 show a comparison of 
experimental and theoretical direct stiffness versus rotor speed for 
the inlet pressure set of table 6 and inlet circumferential velocity 5. 
The figures show that the theory is oversensitive to rotor speed. Some 
of the test data did show an decrease in direct stiffness with rotor 
speed, but not with the sensitivity predicted by the theory. This is 
an improvement over the previous theory of Childs and Scharrer [18] 
which consistently underpredicts direct stiffness. 
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Fig. 79 A comparison of experimental and theoretical direct 
stiffness versus Inlet circumferential velocity ratio at an 
inlet pressure of 3.08 bar and rotor speed of 3000 cpm. See 
table 5 for seal definitions. Teeth-on-rotor (left), teeth* 
on-stator (right). 
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Fig. 81 * comparison of exper mental and theoretical direct 

stiffness versus Inlet circumferential velocity ratio at an 
inlet pressure of 8.25 bar and rotor speed of 3000 cpa. See 
table 5 for seal definitions. Teeth-on-rotor (left), teeth- 
on-stator (right). 
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Fig. 83 A comparison of experimental and theoretical direct stiffness 
versus rotor speed for seal 1 and inlet circumferential 
velocity 5. See table 6 for definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Fig. 84 A comparison of experimental and theoretical direct stiffness 
versus rotor speed for seal 2 and inlet circumferential 
velocity 5. See table 6 for definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Fig. 85 A comparison of experimental and theoretical direct stiffness 
versus rotor speed for seal 3 and Inlet ci referential 
reioclty 5. See table 6 for definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Cross-coupled stiffness 


Figures 86-89 show a comparison of experimental and theoretical 
cross-coupled stiffness versus inlet circumferential velocity ratio for 
the seals defined in table 5. The figures show that, like the test 
results, there is no consistent trend with clearance. The figures also 
show that the theory does an excellent Job of predicting the cross- 
coupled stiffness for the low rotor speed results and a reasonable Job 
for the high rotor speed results. Rotor speed effects are much clearer 
in figures 90-92. Figures 90-92 show a comparison of experimental and 
theoretical cross-coupled stiffness versus rotor speed for the inlet 
pressure set of table 6 and inlet circumferential velocity 5. These 
figures show that the theory predicts reasonably well until higher 
speeds are reached. The theory then predicts a sharp upswing in the 
cross-coupled stiffness at high speeds. This effect is shown by the 
experimental data in figure 90. The larger clearance seals do not show 
this effect at the speeds tested. Perhaps at higher speeds, larger 
clearance seals will show this effect. 
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Fig. 86 A comparison of experimental and theoretical cross<-coupled 
stiffness versus inlet circumferential velocity ratio at an 
Inlet pressure of 3.08 bar and rotor speed of 3000 cpm. See 
table 5 for seal definitions. Teeth-on-rotor (left), teeth- 
* - on^stator (right). 
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Flg . 87 A comparison of experimental and theoretical oroas-coupled 
.stiffness versus Inlet circumferential velocity ratio at an 
Inlet pressure of 3.08 bar and rotor speed of 16000 cpm. See 

-a u@ 5 for seal definitions. Teeth-on-rotor (left), teeth- 
on-stator (right). ' 
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Fig. 88 A comparison of experimental and theoretical cross-coupled 
stiffness versus inlet circumferential velocity ratio at an 
inlet pressure of 8.25 bar and rotor speed of 3000 cpm. See 
table 5 for seal definitions. Teeth-on-rotor (left), teeth- 
on-stator (right). 
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Fig. 89 A comparison of experimental and theoretical cross~coupled 
stiffness versus Inlet circumferential velocity ratio at an 
Inlet pressure of 8.25 bar and rotor speed of 16000 cpm. See 
table 5 for seal definitions. Teeth-on-rotor (left), teeth- 
on-stator (right). 
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Fig. 90 A comparison of experimental and theoretical cross-coupled 
stiffness versus rotor speed for seal 1 and Inlet 
circumferential velocity 5. See table 6 for definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. 91 A comparison of experimental and theoretical cross -coup led 
stiffness versus rotor speed for seal 2 and inlet 
circumferential velocity 5. See table 6 for definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 


157 


xy (N/mm) 


SERL 3 SEAL 3 

SWIRL 5 SWIRL 5 



2000 B000 10000 14000 2000 G000 10000 14000 


4000 8000 12000 1G000 4000 8000 12000 16003 

ROTOR SPEED (CPM) ROTOR SPEED (CPM) 

Fig. 92 A comparison of experimental and theoretical cross-coupled 
stiffness versus rotor speed for seal 3 and inlet 
circumferential velocity 5. See table 6 for definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Direct damping 


Figures 93~96 show a comparison of experimental and theoretical 
direct damping versus inlet circumferential velocity ratio for the 
seals defined In table 5. The figures show that the theory correctly 
predicts an Increase In direct damping for an Increase In clearance for 
the teeth-on-rotor seal. However, the theory Incorrectly predicts the 
same trend for teeth-on-stator seals. This error renders the theory 
suspect when used for teeth-on-stator seals whose geometry differs 
significantly from those tested in this study. However, the theory 
averages an error of 40* for the teeth-on-rotor seals, which Is a great 
improvement over the 75* error of the previous theory of Childs and 
Scharrer [18]. Figures 97— 99 show a comparison of experimental and 
theoretical direct damping versus rotor speed for the inlet pressure 
set of table 6 and inlet circumferential velocity 5. These figures 
show that the theory predicts more speed sensitivity than is shown by 
the experimental data. Perhaps at higher speeds, the test data will 
show the same trends. 

Comparison to theory of [18] 

Figures 100 and 101 provide a brief comparison of the present 
theory to the theory of Childs and Scharrer [18]. Figure 100 shows 
cross-coupled stiffness versus pressure ratio for a teeth-on-rotor 
labyrinth seal at 16000 cpm. The figure shows that the present theory 
follows the experimental data closely while the former theory deviates 
ea pressure ratio is increased. Figure 101 shows direct damping versus 
pressure ratio for a teeth-on-rotor seal at 16000 cpm. The figure 
shows that the present theory follows the experimental data more 
closely than the former theory. 
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Fig. 93 A comparison of experimental and theoretical direct damoln* 
versus Inlet circumferential velocity ratio at an Inlet 
pressure of 3.08 bar and rotor speed of 3000 cpm. See table 5 
deflnltlon3 * Teeth-on-rotor (left), teeth-on-stator 
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Fig. 9* A comparison of experimental and theoretical direct damping 
versus Inlet circumferential velocity ratio at an Inlet 
pressure of 3.08 bar and rotor speed of 16000 cpm. See table 5 
for seal definitions. Teeth-on-rotor (left), teeth-on-stator 
(right). 
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Fig. 95 A comparison of experimental and theoretical direct damping 

versus Inlet circumferential velocity ratio at an inlet 

pressure of 8.25 bar and rotor speed of 3000 cpm. See table 5 

f t° r . w?? al deflnltlons * Teeth-on-rotor (left), teeth-on-stator 
l ri jnt / • 


162 



Cxx (N-sec/m) 




500 

450 

400 

350 

300 

250 

200 

150 


PRESSURE: 8.25 BAR 
SPEED: 16000 (CPM) 



E 

\ 

a 

a) 

w 

I 

Z 


x 

X 

u 


PRESSURE: 8.25 BAR 
SPEED: 18000 (CPM) 



CIRCUM. VEL. RATIO 


CIRCUM. VEL. RATIO 


Fig. 96 A comparison of experimental and theoretical direct damping 
versus Inlet circumferential velocity ratio at an Inlet 
pressure of 8.25 bar and rotor speed of 16000 cpm. See table 5 
for seal definitions. Teeth-on-rotor (left), teeth-on-stator 
(right). 


163 



Cx x (N~ s/m) 


350 


SERL 1 
SWIRL 5 


600 


SERL 1 
SWIRL 5 



2000 6000 10000 14000 2000 6000 10000 14000 


4000 8000 12000 16000 4000 8000 12000 16000 

ROTOR SPEED (CPM) ROTOR SPEED (CPM) 

Fig. 97 A comparison of experimental and theoretical direct damping 
versus rotor speed for seal 1 and Inlet circumferential 
velocity 5. See table 6 for definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Fig. 98 A c caparison of experimental and theoretical direct damping 
versys rotor speed for seal 2 and inlet circumfereSial 

5 * S** table 6 for definitions. Teeth-on-rotor 
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Fig. 99 A comparison of experimental and theoretical direct damping 
versus rotor speed for seal 3 and inlet circumferential 
velocity 5. See table 6 for definitions. Teeth-on-rotor 
(left), teeth-on-stator (right). 
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Fig. 100 A comparison of experimental and theoretical 
results of this report with those of [18] 
for cross-coupled stiffness. 
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Fig. 101 A comparison of experimental and theoretical 
results or this report with those of [18] 
for direct damping. 
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CHAPTER VII 
CONCLUSIONS 

A new analysis utilizing a two control volume model which 
incorporates a solution for the recirculating velocity in the cavity 
and information from a 2-D CFD calculation, has been presented for the 
problem of calculating rotordynamic coefficients for labyrinth seals. 
This analysis was developed to provide both an improved prediction for 
the rotordynamic coefficients and a more detailed model for the flow in 
a labyrinth seal. A seal-test facility has been developed and modified 
for high speed testing for the study of various types of gas seals. A 
method for determining rotordynamic coefficients from experimental data 
has been established, and consistent, repeatable results have been 
obtained. 

A comparison between the CFD results of Rhode and the results of 
new analysis presented in this report support the following 
conclusions: 

(1) The new two-control'- volume model accurately predicts the 
stator wall shear stress for a teeth-on-rotor labyrinth seal cavity. 

(2) The analysis predicts the cavity wall shear stress of a 
teeth-on-rotor seal within 25* of the average of the CFD result. 

(3) The 2-D Jet flow theory used in this analysis accurately 
predicts magnitude of the recirculation velocity along the dividing 
streamline. 

(4) The CFD results show that the mixing length parameter, l, 
used in the equation for the free shear stress is relatively constant, 


for teeth-on-rotor seals, and need not be considered a function of seal 
geometry, as was assumed by Jenny et al. [21]. 

The experimental results of the previous section support the 
following conclusions: 

(1) For teeth-on-rotor seals the direct damping increases as 
clearance Increases; for teeth-on-stator seals, the direct damping 
decreases as clearance increases. 

(2) Direct stiffness and direct damping show little or no 
sensitivity to rotor speed up to 16000 cpm. Cross-coupled stiffness 
shows a sharp upswing at higher rotor speeds, for a teeth-on-rotor 
seal. Cross-coupled stiffness decreases as rotor speed increases, for 
a teeth-on-stator seal. 

(3) Direct stiffness is negative and Increases as clearance 
increases, for both seal configurations. Cross-coupled stiffness 
showed no consistent trend with respect to clearance changes. 

(*•) As clearance decreases, teeth-on-rotor seals become less 
stable and teeth-on-stator seals become more stable. 

The theoretical results of the previous section support the 
following conclusions: 

(1) Theoretical results for leakage underpredict the test 
results presented in this report by about 25). Leakage Increases as 
clearance increases for both seal types. 

(2) Theoretical results for pressure gradient are underestimated 
for teeth-on-rotor seals and overestimated for teeth-on-stator seals. 

(3) The theory correctly predicts that direct stiffness is 
negative and Increases as clearance increases, for both seal 


configurations. The theory Incorrectly predicts an approximately 
quadratic increase In the direct stiffness magnitude (becoming more 
negative) as speed increases. Test results show scant sensitivity. 

(*) The theory accurately predicts an Increase In cross-coupled 
stiffness at high speeds, for a teeth-on-rotor seal. 

(5) For teeth-on-rotor seals, the theory correctly predicts an 
increase In direct damping for an increase in clearance. However, the 
theory incorrectly predicts the same trend for a teeth-on-stator seal. 

(6) The theory Incorrectly predicts an approximately quadratic 
increase in direct damping with running speed. Test results show no 
systematic change in direct damping with running speed. 

(7) A comparison with test results for a teeth-on-rotor seal 
shows that the theory presented In this report does a better job of 
predicting direct damping and cross-coupled stiffness than does the 
theory of Childs and Scharrer [18]. A comparison with the theory of 
Jenny et al. [21] was not possible, as discussed In the THEORETICAL 
DEVELOPMENT chapter of this report. 

In summary, the analysis presented here is considered useful for 
predicting cross-coupled stiffness for both teeth-on-rotor and teeth- 
on-stator labyrinth seals directly. These results were reasonable and 
consistent for the various geometries and operating conditions 
presented. However, the results for the direct damping coefficient for 
teeth-on-stator labyrinth seals were not consistent with test results 
for clearance change effects. This discrepancy renders the analysis 
suspect for predicting these coefficients for teeth-on-stator seals 
whose geometry differs significantly from the seals tested in this 


study. Predictions for the remaining coefficients can be modified using 
the appropriate correction factor from the comparison plots presented 
In this study. 

In the future, if any advances are to be made In the prediction of 
rotordynamlc coefficients for labyrinth gas seals, they will probably 
involve the perturbation of a finite difference solution. The "bulk 
flow" model presented In this report Is too crude to model the complex 
flowfleld present in a labyrinth cavity. 
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APPENDIX A 


GOVERNING EQUATIONS FOR TEETH-ON-STATOR SEAL 


Reduced Equations 


The main difference between the teeth-on-stator equations and the 
teeth-on-rotor equations occurs In the momentum equations. The shear 
stresses acting on control volume I are now the rotor shear stress, \ r , 
and the free shear stress, tj. Similarly, the shear stresses acting on 
control volume II are now the stator shear stress, t 8 , and the free 
shear stress, tj. These difference are evident in the reduce form of 
the continuity and momentum equations given below: 

Continuity I 
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Momentum II 
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where asj and arj are defined as 


(A3) 


as* - (2B+L)/L ; arj -1.0 

The rotor shear stress In the circumferential direction Is now defined 
using the smaller hydraulic diameter and the velocity components of 
control volume I. 
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1 _____ P (Rs |W~W, ) a +U 7* Dh,l\ mr 

fr ■ ~ p/(R9,d)-W, ) Z *U, Z (RSjU-W,) nry J (A 1 )) 

where Dh,j is the hydraulic diameter of C.V. I, defined by 

Dh,i . 2CriL/(Cri*L) (A5) 

Similarly, the atator ahear atreaa in the circumferential direction la 
now defined ualng the larger hydraulic diameter and the velocity 
componenta of control volume II . 

1 //W , a +tV DhjiX® 3 

T s - - p^U7 W* na^ — J (A6) 

where Dh 2 j is the hydraulic diameter of C.V. II, defined by 

Dh al - 2BL/(B-L) (A7) 

The definition of the free ahear stress remains the same. However, 
since CFD results were only available for the teeth-on-rotor 
configuration, the sensitive mixing length ratio l/b may change for a 
teeth-on-stator seal. 

Zeroth-Order Equations 

Continuity: ra ol+1 . i ol (a8) 

Momentum I: ®oi (W,oi“W ,oi-1 ) ■ (tJIo 4 trioari )Li (A9) 

Momentum II: T JioLi - -tsioasiLl (A10) 

Flrat-Order Equations 

The first-order equations remain exactly the same as before. 
Since changes were made in the locations and definitions of the rotor 
and stator shear stress terms, the following changes in the 
coefficients of the first-order equations are necessary: 
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DEFINITION OF THE FIRST ORDER CONTINUITY 
AND MOMENTUM EQUATION COEFFICIENTS 
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APPENDIX C 


SEPARATION OF THE CONTINUITY AND MOMENTUM EQUATIONS 
AND DEFINITION OF THE SYSTEM MATRIX ELEMENTS 
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, X a P 0 i - 

X,«|Pai ♦ w asl 


X| W,8l 

Pci + X a Pci4l * X, W,cl 4 X 7 W aC i - ■o^ici~1 ■ "X,(a4b)/2 

’ x ' [i£ ’ “] - * *■ - *.»]'’* - T57 


>lCl 


*■ X H Psl ♦ X, P S 141 ♦ X, H ,si 4 X 7 W a si - m 0 W ,si-i - 0 


« 
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MOMENTUM II: 


cosiest): Y,„ h * m1 . • |y,« . Y, . ^]p;, 

♦ *0 p cl 4 *• *tcl 4 *. p cl -1 4 *t W»cl - Y,(b-a )/2 

Mn<9 *“ th w *ci - pB? • - pB7 * 

4 ** p ai 4 *• w asl 4 *• p ai -1 4 X? W»al - 0 
Co3(0“iot) : W- Ml • pgii ♦ Y. - Y,.Jp;j . p^2i • M, 

♦ Y% P c j ♦ Y g W aC j ♦ Y # Pcl-1 4 ^7 ^icl m ~Y t (a+b )/2 

aln(e-wt): Y ,u H-, CI - ♦ Y. - Y,^ - ^ w 

4 ^ P 8i 4 ^ ? WjsJ ♦ Y . Pql-1 ♦ Y f Nigl • 0 

Aj~i MATRIX 


aSi 


acl 


*1,2 " A 2 j - A$ t n - Al | # 3 - G, 

A 5 ,2 " A 6 ,1 ’ a 7 ,H “ A 8 ,3 " x s 

a 5,6 ‘ a 6 ,5 “ A 7f 8 “ A8 f7 - -m 0 

a 9,2 • A 10,1 " A 11 f l» - a 12,3 * *• 

The remaining elementa are zero. 

Ai MATRIX 


A 1 t 1 - _A 2,2 - G, at ♦ G a 
A 3,3 " ~ A M ■ G, w ♦ G a 
A 1 ,2 " a 2 , 1 “ A 3 ,*« * A M " °» 

a 5,2 ■ A 6 ,1 - a 7,4 • a 8,3 - X„ 
a 5 , 1 ■ " A 6 ,2 ■ X, w ♦ X a ♦ X,W,i/Ra a 
a 7 ,3 - -a 8 ,H - -X, « ♦ X, ♦ X,W al /Ra a 
a 9 ,2 " A 10,1 - A ii ( n - A 12 ,3 - I* 
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*9,1 * "*10,2 - w ♦ Y, ♦ y a W ai /Rs a 
*11,3 " "*12,H “ "*a <•> 4 T» 4 **W al /Rs a 
*1,5 " ~*2,6 " *3,7 " “All, 8 - C, 



*5,6 - *6,5 - * 7,8 - *8,7 - x . 

*9,6 ■ *10,5 ■ * 11,8 - *12,7 ” *7 
*1 ,9 _ *3,11 “ -* 2,10 “ -AH , 12 “ G* 

*5,9 • "* 6,10 - *7,11 - "* 8,12 " x a p oi /Rs a 
*5,10 * *6,9 * *7,12 ■ *8,11 " x 7 
*9,9 " "*10,10 “ x » “> 4 *i W al /Rs a ♦ 7 a P 0 i/Rs 2 
*11,11 " "*12,12 ■ " x i u ♦ Y, Vf a i/Rs a ♦ 7 t p ol^ Rs * 

* 9,10 ■ *10,9 “ * 11,12 “ * 12,11 - y» 

The remaining elements are zero. 

Aj+T MATRIX 

*1,2 ■ * 2,1 - A 3 f n - Alj t 3 - C 7 
The remaining elements are zero. 
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Fig. D1 Direct stiffness versus rotor speed for seal 1 and inlet 
circumferential velocity 5. See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D2 Direct stiffness versus rotor speed for seal 2 and inlet 
circumferential velocity 5. See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D3 


Direct stiffness versus rotor speed for seal 3 and Inlet 
circumferential velocity 5. See table 6 for pressure raUo 
definitions. Teeth-on-rotor (left), teeth-on-dtator (right). 
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Fig. D5 
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Fig. D6 Cross-coupled stiffness versus rotor speed for seal 3 and 
inlet circumferential velocity 5. See table 6 for pressure 
ratio definitions. Teeth-on-rotor (left), teeth-on-stator 
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circuDfer-enUai velocity 5. See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D8 Direct damping versus rotor speed for seal 2 and Inlet 

circumferential velocity 5. See table 6 for pressure ratio 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D9 
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definitions. T.eth-on-rotor (left), teeth-on-stetor (right)? 
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Fig. DIO A comparison of experimental and theoretical leakage versus 
inlet circumferential velocity ratio for seal 1 at a rotor 
speed of 16000 cpm. See table 6 for pressure definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. Dll A comparison of experimental and theoretical leakage versus 

s^ed veloclty ratio fop 2 at a rotor 

speed of 16000 cpm. See table 6 for pressure definitions 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Tig. D12 A comparison of experimental and theoretical leakage versus 
inlet circumferential velocity ratio for seal 3 at a rotor 
speed of 16000 cpm. See table 6 for pressure definitions. 
Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D14 A comparison of experimental and theoretical pressure 

gradients of seal 1 for a rotor speed of 3000 cpm and inlet 
circumferential velocity 2. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D15 A comparison of experimental and theoretical pressure 
gradients of seal 1 for a rotor speed of 3000 cpm and inlet 
circumferential velocity H. See table 6 for pressure 

definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D16 A comparison of experimental and theoretical pressure 
gradients of seal 1 for a rotor speed of 3000 cpa and Inlet 
circumferential velocity 5. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D17 A comparison of experimental and theoretical pressure 
gradients of seal 2 for a rotor speed of 3000 cpm and inlet 
circumferential velocity 1 . See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D18 A comparison of experimental and theoretical pressure 

gradients of seal 2 for a rotor speed of 3000 cpm and Inlet 
circumferential velocity 2. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 


206 


PRESS. GRRDIENT (BAR) 


9 


SEAL 2, SWIRL 4 
SPEED: 3000 CPM 


SEAL 2, SWIRL 4 
SPEED: 3000 CPM 




• • 



FIXIRL LOCRTION(mm) AXTRL LOCRTION(mm) 


Fig. D19 A comparison of experimental and theoretical pressure 
gradients of seal 2 for a rotor speed of 3000 cpm and inlet 
clrcuDferential velocity 4. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D20 A comparison of experimental and theoretical pressure 

gradients of seal 2 for a rotor speed of 3000 cpm and inlet 
circumferential velocity 5. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D21 A comparison of experimental and theoretical pressure 

gradients of seal 3 for a rotor speed of 3000 cpo and inlet 
circumferential velocity 1. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D22 A comparison of experimental and theoretical pressure 
gradients of seal 3 for a rotor speed of 3000 opm and Inlet 
circumferential velocity 2. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 
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Fig. D23 A comparison of experimental and theoretical pressure 
gradients of seal 3 for a rotor speed of 3000 cpm and inlet 
circumferential velocity See table 6 for pressure 

definitions. Teeth-on-rotor (left), teeth-on-stator (right). 


211 


PRESS. GRADIENT (BAR) 


SEAL 3, SWIRL 5 
SPEED: 3000 CPM 


SEAL 3, SWIRL 5 
SPEED: 3000 CPM 


AXIAL LOCATION (mm) 


AXIAL LOCATION(mm) 


Fig. D2H A comparison of experimental and theoretical pressure 
gradients of seal 3 for a rotor speed of 3000 cpm and inlet 
circumferential velocity 5. See table 6 for pressure 
definitions. Teeth-on-rotor (left), teeth-on-stator (right). 



